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Abstract 

Prior work using gauge/gravity duality has established the existence of a quan- 
tum critical point in the phase diagram of S+l-dimensional gauge theories at finite 
charge density and background magnetic field. The critical theory, obtained by tuning 
the dimensionless charge density to magnetic field ratio, exhibits nontrivial scaling in 
its thermodynamic properties, and an associated nontrivial dynamical critical expo- 
nent. In the present work, we analytically compute low energy correlation functions in 
the background of the charged magnetic brane solution to 4-|-l-dimensional Einstein- 
Maxwell-Chern-Simons theory, which represents the bulk description of the critical 
point. Results are obtained for neutral scalar operators, the stress tensor, and the 
U(l)-current. The theory is found to exhibit a twisted Virasoro algebra, constructed 
from a linear combination of the original stress tensor and chiral U(l)-current. The 
effective speed of light in the IR is renormalized downward for one chirality, but not the 
other, by finite density, a behavior that is consistent with a Luttinger liquid description 
of fermions in the lowest Landau level. The results obtained here do not directly shed 
light on the mechanism driving the phase transition, and we comment on why this is 
so. 



^This work was supported in part by NSF grant PHY-07-57702. 
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1 Introduction and summary of results 



Holography provides a powerful tool for the study of strongly interacting fermion systems, 
including strongly coupled gauge theories. Non-relativistic extensions to systems at finite 
temperature, charge density, and external magnetic field, have applications to physical prob- 
lems in many areas, including in condensed matter, heavy ion collisions, and astrophysics 
[U El [3] . Using holography, thermodynamic functions, correlators, and transport properties 
may be calculated in terms solutions to classical gravity equations. A common setting for 
many of these systems is provided by RG flows in 3+1-dimensional gauge theories, with or 
without supersymmetry. The gravity dual to a wide class of these gauge theories is governed 
by 4-|-l-dimensional Einstein-Maxwell-Chern-Simons theory, whose only free parameter is 
the Chern-Simons coupling k, which is related to the number of chiral fermions in the dual 
gauge theory. 

In the presence of an external magnetic field, the RG flow from 3+1-dimensional gauge 
theory in the UV takes us, in the IR, to an effective 1-1-1- dimensional CFT which is dominated 
by massless fermions in the lowest Landau level. These lowest Landau level fermions are 
only free to move parallel to the magnetic field lines, which leads to the 1+1-dimensional 
character of the low energy theory. This picture is borne out on the gravity side, where the 
RG flow is governed by the purely magnetic brane solution which interpolates between AdSs 
in the UV and AdSa (or BTZ when the temperature is small) in the IR [U [5]. Holographic 
calculations of two-point correlators of the U(l) current and stress tensor further confirm 
the existence of two Virasoro algebras and a single chiral current algebra in the IR [B]. The 
common central charge of the Virasoro algebras equals the Brown-Henneaux value ff\ for 
the AdSs near horizon geometry which, via the Cardy formula, accounts for the specific heat 
coefficiently s/T at low T. Such a specific heat is consistent with a 1+1-dimensional Luttinger 
liquid description [S] of the fermionic theory. The chirality of the current algebra is given by 
the sign of k, while its level is proportional to \k\. Thus, in this regime, a consistent picture 
has emerged in which the behavior of thermodynamic quantities and correlators are found 
to nicely match. 

As the magnetic field is lowered (or equivalently as the charge density is increased), 
the lowest Landau level finds itself increasingly populated. Eventually, we expect qualitative 
changes to occur, both due to occupancy of higher fermionic Landau levels as well as bosonic 
Landau levels, along with any other strong interaction effects at finite density. Numerical 
and analytical holographic calculations, carried out in the dual Einstein-Maxwell-Chern- 
Simons theory, confirm this expectation and reveal a rich phase structure as a function of 
the Chern-Simons coupling k [9l [TOl [11] (see [121 US] for other examples of holographic phase 



^In the boundary CFT only dimensionless combinations have an invariant physical meaning. Dividing 
by a suitable power of the charge density p, the dimensionless ratios for the magnetic field B, the entropy 
density s, and the temperature T will respectively be denoted hy B = B/ p^^^, s = s/ p, and T — T j p^l"^ . 
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transitions involving magnetic fields and finite density.) Specifically, for k > 1/2, the low 
T specific heat coefficient was derived exactly as a function of the magnetic field B, and is 
given by. 



s 
f 



TT 



6 53 



53 



B > Br 



XI) 



While the constancy of s/T as T — )■ is still as in Luttinger liquid theory, its dependence on 
B is not. As the magnetic field approaches the critical value B^ (which depends only on k), 
a quantum critical point was discovered at B = Be- The scaling behavior of s at i? = i?c was 
obtained numerically in [10], and subsequently established by analytic calculations [TT], both 
for > 3/4. As reviewed in Appendix |A| the picture was completed since by a numerical 
calculation of the scaling in the range 1/2 < k < 3/4. Assembling these data, we have. 



a 



1 

3 
k 



l<k<l 
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This result signals further strong deviations from canonical Luttinger liquid behavior. The 
full scaling function, accounting for variations in both T and 13 in the critical region, was 
also computed analytically in [TT] for the range k > 1, and agrees with the numerical results 
of jlOj. 



Remarkably, as discussed in [TO], the critical exponent 1/3 arises also in the Hertz/Millis 
theory [lH [15], [16] of metamagnetic quantum phase transitions, in which a variation in the 
magnetic field across a non-zero critical value produces non-analytic behavior in the specific 
heat, but no change in symmetry. Metamagnetic quantum criticality has been observed 
[T7] in a wealth of recent experiments on Strontium Ruthenates, where the entropy density 
has been mapped out through an extensive region of the temperature and magnetic field 
parameter space. Interestingly, as the critical point is approached the specific heat coefficient 
is found to diverge linearly in 1/{B — B^), just as in (1.1). 

The purpose of the present paper is dual. On the one hand, we wish to extend to the 
case of finite charge density the holographic calculations of various correlators which have 
been evaluated in [S| for vanishing charge density, and understand the relation with the low 
T thermodynamics already established for finite charge density. On the other hand, if the 
holographic calculations of thermodynamics are really relevant to metamagnetic quantum 
phase transitions, then obtaining the corresponding holographic correlators should allow one 
to address transport problems in metamagnetic materials close to quantum criticality in a 
model with only a single free parameter k. The results of this paper will constitute a modest 
first step in both directions, and will be summarized below. 
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1.1 Summary of results 

In this paper, 2-pomt correlators will be evaluated in the IR regime for a neutral scalar, for 
the U(l) current, and for the stress tensor, in the presence of the charged magnetic brane 
solution. The method of overlapping expansions may be justified for small momenta, and 
will be used throughout. 

Calculating the 2-point correlators of a neutral scalar in the IR regime, we identify 
effective left- and right-moving excitations with different propagation speeds. While left- 
movers propagate at the speed of light (with respect to the standard boundary metric) as 
they do for the purely magnetic brane, the speed of the right-movers is less than the speed of 
light in the presence of non-zero charge density. We compare these holographic predictions 
for the dynamics of chiral excitations with those of Luttinger liquid theory [8j , in which the 
speed is also predicted to be renormalized downward. Finally, we argue in favor of a relation 



between the entropy density scaling of (1.2) for 1/2 < k < 3/4 and the dynamical scaling 
properties of the scalar two-point functions. We leave unexplained, however, the precise 
nature of the crossover to the critical exponent 1/3, and the precise value of the coupling k 
at which this cross-over should occur. 

Calculating the 2-point functions of the U(l) current J'± and the stress tensor T±± reveals 
that their IR behavior is governed by two Virasoro algebras, one of which is twisted with 
the chiral current algebra. The positive chirality correlators are given by the following 
expressions, valid for all A; > 1/2, 



kc 1 



27r2 {x+ - y+Y 
kc^ 1 



27r2 (x+ - y+f 

n4x^)n4y'-)) ~ ^ ^ +^2 + ^ r / +^4 (1-3) 

' 2tt^ (x+ — y+)'' 8tt^ (x+ — y+j* 

where c is the Brown-Henneaux central charge, and fi is the chemical potential. The non- 
vanishing negative chirality correlator is given by, 

/ STT^ [x — y ) 

These correlators have higher order IR corrections, whose quantitative evaluation would 
require going beyond the overlapping expansion method. Qualitatively, however, we know 
that some of these higher order corrections arise from double and higher trace operators 
which are generated by the RG flow towards the IR [IHl 112], just as they were shown to arise 
in the purely magnetic case in j6]. Upon forming the combination 

ri^ = r++ + (1.5) 
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the chiral current algebra of J7V decouples from the Virasoro algebra of 7^+, and we are left 
with two decoupled Virasoro algebras, and one decoupled chiral current algebra, just as we 
had in the case of the purely magnetic case. These results are valid for all k > 1/2, and 
constitute perhaps the most important result of this paper. The net effect of turning on a 
non-zero charge density (or equivalently, a non-zero chemical potential fi) is simply to twist 
the positive chirality Virasoro algebra by the admixture of /ij7+, as well as renormalizing the 
speed of propagation. These results confirm the presence of a unique underlying CFT for 
all values of the magnetic field, but offer no immediate insight into the dynamical scaling 
properties of the entropy density in (1.2). 



1.2 Organization 

The remainder of this paper is organized as follows. In section |2j we review the Einstein- 
Maxwell-Chern-Simons theory and the salient properties of its charged magnetic brane solu- 
tions. In section [3| we calculate the 2-point correlators of a neutral scalar field of arbitrary 
mass, and discuss some of their immediate physical consequences in section |4} In section [5} 
we outline the general structure of the 2-point correlators of the U(l) current and the stress 
tensor. In section |6] we set up the technical calculation of the linear fluctuations in the gauge 
field and metric using overlapping expansions in light-cone gauge for the near region to LMN 
and Fefferman- Graham gauge in the far region. The detailed calculations of the twisted Vi- 
rasoro terms in the current and stress tensor 2-point functions are derived in sections [7] and 
[8j A brief discussion of some open questions is presented in section [9j In Appendix |X] we 
review the behavior of the low temperature specific heat exponent as a function of k. In 
Appendix [B| perturbative and WKB approximations to the scalar correlators are developed. 
In Appendix [C} the full linearized fluctuation equations around the charged magnetic brane 
solution are derived in light-cone gauge, while the gauge transformations required in the 
overlap region are constructed in Appendix [D] 
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2 Review of the charged magnetic brane solution 



The action of five-dimensional Einstein-Maxwell-Chern-Simons theory with a negative cos- 
mological constant and Chern-Simons couphng k was given in [6]. Here, we will need the 
corresponding Bianchi identity dF = and field equations, 

= d^F + kFAF 
Rmn = ^Qmn + T^F^^ FpggMN — '^FmpF^^ (2.1) 

Throughout, we will take A; > 0, since a sign reversal of k is equivalent to a parity trans- 
formation. For the supersymmetric value k = 2/a/3, this system is a consistent truncation 
for all supersymmetric compactifications of Type IIB or M-theory to AdSs (see [201 



Solutions of (2.1) are guaranteed to be solutions of the full 10 or 11 dimensional supergravity 
field equations, and to provide holographic duals to the infinite class of supersymmetric field 
theories dual to these more general supersymmetric AdSs compactifications. 

2.1 Boundary stress tensor and current 

We will be considering asymptotically AdSs solutions, for which the metric and gauge field 
admit a Fefferman- Graham expansion [23j. Introducing a radial coordinate p, defined such 
that the AdSs boundary is located at p = oo, the metric takes the asymptotic form 



ds^ = — - + gny{p,x)dx^dx'' 
4p^ 



g^P^x) = pg^Kx)+g^^{x) + -gf}{x) + ^^g^^{x) + --- (2.2) 



P p 



while the gauge field takes the asymptotic form. 



A = A^{p,x)dx- 
A,ip,x) = A^;\x) + U^;\x) + --- (2.3) 

Here x^ = {p,x^), with p = +, — , 1,2. The coefficients g^^J, g^^^\ and the trace of g^^J are 
fixed by the Einstein equations to be local functionals of the conformal boundary metric gj^J. 

The boundary stress tensor and current [2U [25] are defined in terms of the variation 
of the on-shell action with respect to gjj^J and A^^^ respectively (see [6J). In terms of the 
Fefferman-Graham data the result is, 

4nG,T^U^) = ^?(tna;) + local 

2'kG^J^{x) = A(f)(x) + local (2.4) 
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Indices are raised and lowered using the conformal boundary metric gj^J. The local terms 
denote tensors constructed locally from gj^J and A^^\ In this paper, as in [6], we are interested 
in computing two-point correlation functions of operators at non-coincident points, to which 



the local terms in (2.4) do not contribute. Henceforth, we will drop all local terms. 



2.2 Ansatz, Reduced field equations, and symmetries 

The charged magnetic brane solution was constructed in [TT]. It belongs to a class of so- 
lutions characterized by a constant background magnetic field, rotation invariance around 
this magnetic field, and translation invariance in all boundary directions x^. The Ansatz is 
conveniently expressed in the so-called "LMN" gauge, 

A = bx^dx^ + A+dx^ + A_dx~ (2.5) 

The index i = 1,2 labels the coordinates of the plane orthogonal to the magnetic field. The 
magnetic field strength b = F12 is constant, and by rescaling x^'"^, may be set to a convenient 
value, b = ^/S. The functions L, M, N, V, A± depend only upon the holographic coordinate r. 
The reduced field equations were given in [T^, with the identifications E = A'^, P = —A'_, 

Ml (^{NE + LP)e^^y + 2kbP = (2.6) 

M2 ({LE + MP)e^^)' -2kbE = Q 

El L" + 2V'L' + A{V" + 1/'^)L - APE = 

E2 M" + 2V'M' + A{V" + V''^)M + 4^^ = 

E3 A^" + 2V'N' + A{V" + V''^)N + 4p2 = 

E4 f{V'f + fV - 6 + ^(L')^ - ^M'N' + b'^e-^^ + MP^ + 2LEP + NE"^ = 

fV /" + Af'V' + 2fV" + Af{V'f - 24 = 

Throughout, primes will denote derivatives with respect to the holographic coordinate, and 
we will use the abbreviation, 

f = L^ - MN (2.7) 



Notice that LMN gauge of (2.5) differs from Feffer man- Graham gauge of (2.3) only in the 
structure of the term proportional to dr'^ . Finally, under constant SL{2, R) transformations 
of the coordinates x^, 



A-M;_ ) AeSL{2,R) {21 
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the Ansatz (2.5) is invariant provided the fields transform as follows, 



while the field V and the combination / are invariant. 

2.3 Charged magnetic brane solution 

The basic charged magnetic brane solution satisfies the consistent restriction, 

N = P = (2.10) 

As a result, the functions L = Lq and V = Vq satisfy the reduced field equations El, E4, 
and fV of the purely magnetic background solution [6J. The E4 equation is a constraint 
whose derivative vanishes by virtue of equations El and fV. Although no analytical solution 
is known, the relevant parameters governing the small and large r asymptotics are available 
numerically, and this will suffice for our purpose. The remaining functions = Aq and 
M = Mq were solved for in [11], and are given by, 

kb 

M„H ^ Ur){-^^-,t,££^) (2.11) 
where ip is defined by, 

, , r dr' , , 

The parameters cv,ce, and aoo entering these functions specify the asymptotic behavior of 
these functions, as we will now spell out. 

2.4 Asymptotic AdSs behavior 

As r — )■ oo, the solution approaches AdSs, and we have the following asymptotics. 

Loir) = 2{r-ro)-^ + Oir~') 
g2v;,(r) ^ cy(r-ro) + 0(r-i) 

. , . CyCE Ce _2\ 

kb r 
Mo(r) = -^r + C(r°) 

^(r) = --^ + 0(r-2) (2.13) 
ley r 
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The interpretation of these constants was given in [T^ . The constants cy and ro characterize 
the purely magnetic solution. The asymptotics EQ{r) ~ ce/t"^ of the electric field Eq = 
shows that ce is proportional to the electric charge density p of the dual state, while the 
asymptotics of Aq gives its chemical potential fj,. The precise relations are as follows, 



1 



where the coefficient a^o arises from the asymptotics of the metric function M( 



0- 

2.5 Asymptotic near-horizon behavior 

As r — 0, the solution approaches a family of exact solutions for which V is constant, 

F = bdx^ A dx"^ + q r'^'^dr A dx~^ 

dr"^ ( 2(7^2'= \ 

ds^ = 7717 - aor + — ^ ^ (dx^f + 4brdx^dx~ + dx'dx' (2.15) 

4o^ Y k[2k — 1) / 

The corresponding functions are given by. 

Loir) = 2br + 0{r'+n ^ = _ J + ^ 

^2Vo{r) ^ 1 + 2r" + C(r2'^) 
Ao{r) = Ir' + Oir'^^n 

rh 

Mo(r) = -aor + C(r^+") 
V'(r) = ^ + ^-0 + 0(0 (2.16) 

The constant ipo is a property of the purely magnetic solution only. On the other hand, q 
and ao are associated with the presence of non-zero charge density. 

The metric appearing in the {r,x~^,x~) directions is the same as a 2+1-dimensional 
Schrodinger solution [2ni EZ] , or alternatively a null- warped solution 



2.6 Relations between asymptotic parameters and regularity 

The key asymptotic data in the large and small r asymptotics are related as follows, 

_ CyCE ofcfeVo 

q ^ 

«o = «^-1644J(fc) '^^^^ = 2kl Lo(r)=^e2^°W ^^'^'^^ 
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They link the large r and small r behavior of Aq and Mq. The integral J{k) is convergent 
for k > 1/2, and is positive, so that we must have < Ooo as long as g 7^ 0. 

Regularity of the solution requires Mo(r) < for all r > O.footnoteMore precisely, if 
this condition is not obeyed then the zero temperature solutions studied here cannot be 
obtained as limits of smooth, finite temperature solutions. When g 7^ 0, the near-horizon 



metric (2.15) is regular only when k > 1/2, a condition we will assume throughout. Since 



Mo(r)/Lo{r) + aoo/{2b) vanishes at r = 00, and is monotonically increasing for decreasing 



r, regularity also requires < Oq. The relation between a^o and ao was shown in [TT|^to be 
governed by the (normalized) magnetic field B = B/ p^l"^ , where B = 2b /cy, 

^ = 1-^ Bl^^im (2.18) 

Thus, regularity of the charged magnetic brane solution forces B > B^, with the critical 
limit B = Be corresponding to a quantum phase transition [TT]. For B < Be a new branch 
of solutions opens up; these were studied numerically in [TT], but no analytical results are 
currently available. 

2.7 The action of SL{2, R) 



SL{2, R) transformations which preserve the consistent restriction P = = of (2.10) may 
be parametrized as follows. 

They will map charged magnetic brane solutions with different parameters into one another. 



Denoting the transformed quantities by tildes, as in (2.9), we find 



q = Xoq 
tto = Aq(q;o + 46A) 

aoo = Xl{a^+4bX) (2.20) 

We see that the parameter Aq rescales the charge q, while the parameter A changes the 
normalized magnetic field. In particular, all charged magnetic brane solutions may be ob- 
tained as SL{2, R) maps of the critical solution for which = 0, B = B^, and a^o = ac = 
16cyc|;J(fc), the magnetic field being given by, 

= ^ (2.21) 



^The coefficients aoo and ao were denoted respectively by a and a in [IT] . 
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Turning on a non-zero temperature, as we did in [TT], or non-zero momenta, as we will 
do here, will make these solutions physically inequivalent, as the SL{2, R) transformation 
will now also change the temperature and/or the momenta. The practical advantage is 
that, without loss of generality, we can restrict to studying the solution corresponding to 
the critical magnetic field B = Be, since the solution for other (finite) values of B may be 
derived from the critical solution by applying the above 5'L(2,R) transformation. Thus, 
henceforth, we set = 0, and ctoo = Oc- 
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3 Correlators of a neutral scalar field 



In the case of the purely magnetic background solution, the fluctuations of the gauge fleld 
decoupled from the metric fluctuations to linearized order, so that the correlator between the 
f/(l) current and the stress tensor vanished [6]. For the charged magnetic brane solution, 
no such decoupling takes place, a fact that renders the problem more involved. A simpler 
warm-up case is provided, however, by studying the fluctuations of a free neutral scalar fleld 
in the presence of the charged magnetic brane solution. Although the neutral scalar will only 
probe the metric of the background solution, and not its gauge fleld, its study will provide 
valuable information for scalar correlators that would be much more difficult to obtain for 
the full stress tensor and current correlators. 



3.1 Reduced Field Equation 



We consider a free neutral scalar fleld $ with mass m. It obeys the wave equation, 

1 



'M 







(3.i: 



The background being invariant under translations in x^, and equation (3.1) being linear, 
it will suffice to consider plane wave solutions with momenta = (p+,p_,pi,p2) 
will be interested only in the x^-dependence, we set pi = P2 = 0, which leaves. 



Since we 



[r,P±)e 



ipx 



(3.2) 
P+P-- 

In the background of the charged magnetic brane solution, the wave equation ( |3.1[ ) reduces 
to, 



where we continue to use the notations of [6j, namely, px = + p^x , and 



2p2 



+ 



plM 
L2 



m 







(3.3) 



Throughout, we will often use the abbreviation e^^ = K. Since the metric is unperturbed in 
this system, we have dropped the subscript on the functions L, M, K. No analytic solution 



is known to the reduced scalar equation (3.3), which comes as no surprise since the functions 



K, L, M themselves are not known analytically. Nonetheless, for small momenta (precise 
conditions will be spelled out later), we may use the method of overlapping expansions 
which was used already in |6] to calculate the entropy density of the charged magneti c br ane 
solution for low temperature, as a function of B. Following this method, we solve (3.3) in 



the low momentum limit, by matching the solutions in overlapping near and far regions, 
corresponding to small and large r respectively. As announced in section 2/7, we restrict 
to analyzing the critical solution corresponding to B = Be, since all other solutions may be 



recovered by the SL{2, R) transformations of section 2.7 
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3.2 Solutions in the far region 



In the far region, momentum dependence in (3.3) is to be neglected, which requires <^ r 
and p"^ r. We want the far region to extend down to r <^ 1 so that the far region will 
overlap non-trivially with the near region; this requires <^ 1 and <^ 1. For k < 1, 
there is, however, an additional requirement: since M(r) ~ r'^^ for r <C 1, the above overlap 
condition will hold provided we have p2 ^,2^-2 ^ x as well (this condition is automatic for 
> 1 in view of the requirement <^ 1 already imposed earlier). Note that it is possible 
to satisfy all these conditions by making p_ small, and yet allowing p+ to be of order 1. 

Assuming that the above conditions are satisfied, all momentum dependence in (3.3) may 
be dropped, and we are left with the reduced equation in the far region. 



LV + ^(irL2)V-mV 
K 







(3.4) 



which involves only the data of the purely magnetic background solution. Clearly, no ana- 
lytical solution for is available, since K and L themselves are known only numerically. The 
data we need from the far region, however, depend on only a small number of parameters 
which may be computed numerically, if so desired. 

To see how this works, consider (3.4) for r ^ 1, where L{r) ~ 2r and K{r) ~ Cyr, 







For r ^ 1, the general solution is given by, 

(j)(r) = a^r'^^ + a-r^~ iz-t = — 1 ± \Jl + m? /A 

For r ^ 1, but still in the far region where momenta are dropped, we have L(r) 
K{r-) ~ 1, so that (3.4) becomes, 

12r V + 24r0' - = 



The general solution of (3.7) is 

Mr) = b>r^^~^ + 6_r" 



1 + mV3 



(3.5) 

(3.6) 
2br and 

(3.7) 
(3.8) 



By linearity of (3.4), there exists a linear relation between the coefficients a± and b±, which 
we may express as follows. 



0+ 
a_ 



S 



6_ 



S 



The correlator of the scalar field $ is given by the ratio, 

^„ 5_+ + S__(6_/6+) 



G 



«+ S++ + S+-{h_/h+) 



(3.9) 



(3.10) 
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The matrix S depends only on the data of the purely magnetic solution, and on the mass m 
of the scalar, but not on any momenta. The ratio b^/b^ will be governed by the dynamics 
in the near region, and will depend on the momenta p±. 



3.3 Solutions in the near region: purely magnetic case 

We begin by solving for the behavior of the scalar in the purely magnetic case, since this may 
be done analytically and simply. The near region is defined by r <^ 1, so that L{r) ~ 2br 



and K(r) ~ 1. Equation (3.3) then takes the form 

-.2 



12r>" + 24r(j)' 



— — m^cj) = 
or 



By the change of variables r 



and the redefinition of the function 



equation (3.11) is transformed into a modified Bessel equation, 

Z (p + Z(p 







(3.ii: 



Z0[Z) 



(3.12) 



Of its two linearly independent solutions, Ky vanishes at the horizon z = oo, while I^, diverges 
there. Retaining only the regular solution, and expressing it in terms of 0(r), we obtain. 



(r) 



2 sin z/TT 
i^Ki, 




(3.13) 



In the r ^ part of the near region, is given by the following approximate behavior. 



463 



1 



r(l-z/) \Ab^) r(l + z/) \Ab^^ 

This formula was derived using the familiar definition. 



1 

7" 2 2 



KJz) 



2 sminu] 



{I^,{z) - h{z)) 



(3.14) 



(3.15) 



as well as the 2; — )■ asymptotics of Iu{z), given in (B.7) of Appendix B. 



3.3.1 IR behavior of 2-point correlator: purely magnetic case 



The far and near regions overlap when r <^ 1, so that the solution (3.8) valid in the 



far region must match the solution (3.14) in the near region. We see that their functional 



dependence on r indeed coincides, which allows us to extract the key ratio, 

b_ _ T{l-u) /p^y 
b^ ~ ~r(i + i^) 1,463 J 



(3.16) 
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Assuming a fixed value of v away from the positive integers, and that ^ 1, we find 
6-/6+ ^ 1. The scalar two-point function of (3.10) may then be expanded in powers of 

6-/ 



and the leading non-trivial momentum dependence identified, 



G{V) 



rfi-z/) 



rri 



22,. ^3!. 



(3.17) 



The first term in G(p) produces a delta function in position space, which is local, and may be 
dropped. The determinant factor det(S') is proportional to the Wronskian of the differential 
equation (3.4), evaluates to det(S') = "ivicy^'^iy^ — vS)~^^ and thus never vanishes. As a 
result, the leading IR behavior of the scalar Green function is given by . 



3.4 Solutions in the near region: charged magnetic brane case 

In the near region r <^ 1 we have iiT ~ 1, L ~ 26r, and M is given by 



M(r) = -auT 



2k 



2g2 



(3.18) 



k{2k - 1) 

where we are restricting to the critical solution and set ao = 0. All momentum dependence 



in (3.3) now needs to be retained, resulting in the equation 

^2 



12r>" + 24r0' - 



P , P-ttM 2k-2, 2 

-0-— r 







Changing variables to r = 1/2;^, and redefining 0(r) = z(j){z), transforms (3.19) into, 

-.2 



z^6" + z(j)' 



P^ 2 , 2 I P-'^M 4_4fc 

—z + z/ H z 

36 







(3.19) 



(3.20) 



Here, u is given by (3.8). For q = um = 0, this equation coincides with (3.12), which has 
already been solved earlier. For p+ = 0, but p_ 7^ 0, equation (3.20) may be transformed 
into a modified Bessel equation, and solved analytically as well, for all k. 

Henceforth, we will assume q,p± 7^ 0. Rescaling z and the scalar field as follows. 



P 



[z) = ip{x) 



transforms (3.20) into. 



x'^if' 



" + xip'-{x^ + iy^ + ex'-"')ip = 
17 



(3.21) 



(3.22) 



Thus, for fixed k and z/, and up to the above rescahng of z, the scalar field equation in the 
near region intrinsically depends only on a single combination ^ of the momenta, given by. 



)l~3k 



\2k{2k 



(3.23) 



Let (p{x) denote the solution to (3.22) that is smooth at the horizon, x — )■ oo. In the overlap 
region, we have ^ r ^ 1 which translates to p ^ a; ^ 1 in terms of the coordinate x, 
and the function ip{x) takes the simplified form, 

^{x) ~ c^O^" + c+iO^-" (3.24) 

The dependence on k and u of coefficients c±{^) has not been exhibited here. Transforming 
back to the holographic coordinate r, we have. 



0(r) ~ c_(0 



P 

63/2 



+ ^+(e)(^) ri-^ (3.25) 



3.4.1 IR behavior of 2-point correlator: general form 



In the overlap region, where we have p <^ r <ti 1, the functional form (3.25) of 0(r) in 
the near region coincides with the functional form (3.8) in the far region. Matching the 
corresponding coefficients gives. 



(3.26) 



Again, since <C 1 we have h-/h+ ^ 1, so that the two-point correlator for the scalar may 
be reliably approximated as follows, 

S_+ det(^) 



G{p) 



The ^ — )■ limit reduces to the purely magnetic case so that we have Cj,(0) 



(3.27) 



C,, with Cu 



given in (3.17). Also, the p+ — )■ limit, while keeping p_ fixed, should be smooth. For A; > 1, 



this limit implies ^ — ?■ 0, and again leads to the purely magnetic case. For /c < 1, the p+ — 
limit implies ^ — )■ oo, and smoothness of G then requires, 

CAi) ~ as ^ ^ oo (3.28) 

The precise coefficient can be computed analytically but will not be needed here. The 
remaining p_ dependence of 6-/6+ is then of the form {j)JyI'^^^^\ as may also be seen 



directly from scaling arguments on the original equation (3.20) with p 



P^ 



0. 



Beyond these asymptotic limits, it does not appear possible to solve equation (3.22) for 



the near region behavior of the scalar field analytically for generic value of k. Besides the 
trivial k = \ case, there is one special value A; = 3/4 where a complete solution may be 
obtained in terms of Whittaker functions. 
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3.5 The solvable special case of /c = 3/4 



For the special value A; = 3/4, equation (3.22) is related to the Whittaker equation [22] • The 
solution to (3.22) which vanishes for large x is given by, 

e ^ 



1 



u+ ^,l + 2u, 2x 



where U is the confluent hypergeometric function in the notation of 
region, its behavior simplifies to. 



x'^ 



X 



V{l + 2u)V{\ 



(3.29) 
In the overlap 

(3.30) 



z/ + f) r(i-2z/)r(i + z/+f) 

where we have omitted an immaterial multiplicative constant common to both terms on the 



right hand side of (3.30). In terms of the holographic coordinate r, we find. 



r(i + 2z/)r( 



r(i 



2u)V{\ 



As a result, we find the following formula, 

^ = a(Op'^ a(o 



r(i-2z/)r i + i/ 



r(i 



2u)v{\ 



(3.31) 



(3.32) 



One may check that the ^ — ?■ limit of the above expression reproduces the general result 
Cj,(0) = Cy of (3.17), by using the duplication formula for the F-function. 

There are three more special cases, k = 5/8, k = 7/8, and k = 3/2 where the solutions 
to the differential equation (3.22) are expressible in terms of Heun functions. As of yet, we 
have found insufficient information on the asymptotic behavior of these functions to make 
use of them in the generation of the scalar two-point correlator for these values of k. 



3.6 Comparing expansions with the exact /c = 3/4 solution 

Beyond obtaining the exact near-horizon solution for A; = 3/4 it is possible to derive the 
functional form of the expansion in powers of ^ of C,y{C,), through a combination of perturba- 
tive and WKB methods. These calculations are carried out in Appendix [B| and their results 
may be summarized as follows. The perturbative expansion in ^ is valid when ^ <^ 1, while 
the WKB expansion holds when /i -C 1, 

- < < 1 in powers of ~ 

2 p_^_ 

1 l/(2fc-l) 

- < k < 1 in powers of ^ ~ ( -^jr- 1 (3.33) 
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Both expansions may be applied to the = 3/4 special case, where we have an exact solution. 
For k = 3/4, the above expansion parameters become. 



f2 



3/2 
1/2 



1 



(3.34) 



The exact solution of (3.32) for A; = 3/4 



r(i + e- 



(3.35) 



clearly admits a Taylor series expansion in powers of for small ^. It admits an expansion 
for large ^, which may be computed using the following series expansion of the F-function, 



0{y 



-2> 



1 e 

- + — 

2 2 



(3.36) 



showing that the expansion parameter is h ~ 1/^^ for large ^. This result matches the 



prediction of WKB given in (3.33), and (3.34). 



Note that, contrarily to the case of the current correlators in the purely magnetic brane 
solution, the corrections in powers of ^, either for large ^ or for small ^, do not readily admit 
an interpretation in terms of double trace deformations of the action and operators, because 
the expansions in powers of ^ or 1/^ are not given by geometric series. It would be interesting 
to investigate whether the structure seen here can be accounted for by multi-trace operator 
perturbations, or requires a novel paradigm. 

We close with a discussion of the case k > 1. The combination now vanishes uniformly 
as p± — )■ 0. It is multiplied by a factor which becomes singular at 2 = when k > 1. 

This singularity is beyond the domain of validity of the near-horizon region, in which we are 



considering equation (3.22). Thus, it is more suitable to start from equation (3.19), with 



r ^ 1 and to carry out a power expansion in p_aM- The lowest order contribution to the 
series is just the purely magnetic solution. 
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4 Physical Applications 



In this section we discuss two immediate physical apphcations of the calculations of scalar 
correlators, given in the previous section. A first concerns the connection between the critical 
scaling law of the entropy density of (1.2). A second deals with the velocities of the chiral 
excitation modes, their dependence on the magnetic field, and similar effects encountered in 
Luttinger liquid theory. 



4.1 Dynamical scaling exponent and entropy scaling form 

We begin by commenting on the dynamical scaling exponent that emerges in the near- 
horizon geometry, and its relation with the thermodynamic scaling relation (1.2) for the 
entropy density. The metric of the near-horizon geometry at the critical point is of the 
Schrodinger type and given by, 

ds^ = — - — Adx+f + Abrdx+dx~ + dx'dx' (4.1) 

Ar'^ k{2k — 1) 

As shown already in [TT], it is invariant under the following scale transformations, 

r Ar a;+ A~''a;+ x' X'^'^x' (4.2) 

This dynamical scaling symmetry is reflected in the IR behavior of the correlators of the 
scalar field through the dependence on only the combination, 

e ~ ip-)\p+)'-' (4.3) 



which is invariant under the scalings of (4.2). Reverting to coordinates x"*" — )■ t and x~ — )■ X3, 
and frequency/momentum p+ — w and p- — ?■ k, used in pTl] to study thermodynamics, we 
readily identify the standard dynamical scaling exponent, defined by w ~ k^, as 

k 



k 



(4.4) 



Scaling arguments alone then predict the scaling relation for the entropy density s as a 
function of T in ci space-dimensions, s ~ T"^/^. Assuming that the effective IR theory is 
indeed a 1+1-dimensional CFT, as well as the value of the dynamical scaling z of (4.4) 
would lead to 

s ~ f " a= (4.5) 

k 



Comparison with (1.2) reveals agreement with the full gravity calculation only for 1/2 < 
k < 3/4, but disagreement for k > 3/4. Thus, the key problem in explaining (1.2) becomes 
to understand the transition across k = 3/4. The near-horizon expansion in powers of ^ will 
also result when we will study the 2-point correlators of the current and stress tensor. 
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4.2 Speed of chiral excitations 

The scalar correlator computed in the background of the charged magnetic brane may be 
used to calculate the velocities at which the low momentum scalar excitations propagate in 
the dual CFT. For k > 1, the low momentum limit \p±\ <^ 1 uniformly corresponds to ^ <^ 1. 
For 1/2 < k < 1 the limit is not uniform, however, and we need to distinguish two regimes, 
loosely characterized by <^ and ^ \p+\, and more precisely specified by ^ -C 1 
and ^ ^ 1 respectively. Up to local terms, and up to an overall constant, the Green function 
of (3.27) in both regimes is given as follows. 



l-fc , , k 



e>l G(p) ~ (p_)i^ (^l + ci(p+)^(p_)~^ + 0(r)j (4.6) 

Since the second line arises only in the interval l/2</i;<l, itis instructive to specialize to 
the solvable value k = 3/4, where G{p) is given by, 

G{p) ~ (p.r + c,ip+)kp-r~'^ +■■■ (4.7) 

Taking the Fourier transform gives the propagator as follows. For 1/2 < k < 1 and small 
compared to well as for k > 1 and both large, we have. 

For 1/2 < k < 1 and x~^ large compared to x~ , we have, 

Gix^) 3^.1 (4.9) 

(x+)2(x-)^''-2 



where the leading (p-)^*^ term in (4.7) is local and may be dropped. 

As is manifest from the factorized form of these contributions, the propagation of left- 
movers and right-movers are independent from one another. To calculate the velocity of 
propagation of each chiral mode, we need to focus on the normalization of the boundary 
metric which, for the critical solution B = Be, is given by, 

ds'^ = 4:dx+dx- +I3(dx+f + dx'dx' (3 = -^ (4.10) 

b 



Thus we see that the propagators of (4.8) and (4.9) were in fact expressed in an oblique 
coordinate system. To find the physical velocities of propagation, we need to convert to a 
set of standard coordinates, x^, defined by 

x~^ = x"*" x~ = x~ + /3x'^ (4.11) 
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in which the boundary metric is given by, 



ds^ = Adx^dx + dx^dx'' 



(4.12) 



Introducing cartesian coordinates t, X3 by x^ = {x^ ±t)/2, the speed of hght is set to 1. The 
propagation of the chiral modes follows the singularities of the propagator. For the mode 
x~^ = 0, we have x~^ = x^ + t = 0, corresponding to a left-moving chiral excitation propagates 
at the speed of light. For the mode x~ = 0, we have x~ — (3x~^ = 0, or (1 — f3)x3 = (1 + f3)t, 
which corresponds to a right-moving chiral excitation propagating at the speed 



1 + 
1-/3 



(4.13) 



Since /5 < 0, the right-mover propagate slower than light. We may obtain a formula for the 
velocity in terms of intrinsic scale invariant physical observables i n the boun dary CFT by 
expressing a^o as follows, a^o = ic^CyB^ /3, using the formulas of (2.17) and (2.18). Using 
next the expression for the chemical potential /i and its scale invariant form fi = jj,/p^^^, we 
further express cyCg in terms of /i, resulting in aoo = Sbk"^ fi"^ B'^ , and the velocity. 



1 - 2k'fx'E 
1 + 2k^p,^B 



(4.14) 



4.3 Comparison with Luttinger liquid theory 

Luttinger liquid theory extends the Fermi liquid picture in 1+1 dimensions by including all 
(marginal) four-Fermi interactions. We begin by reviewing the Luttinger liquid construction. 
As a finite charge density is turned on, a Fermi surface develops, which in 1+1 dimensions 
consists of just two points with opposite momenta ikp, as depicted in Figure 1. 



E 





















Figure 1: Fermi "surface" and 1+1-dimensional chiral fermions V'l and V'r 

At those points, we have chiral fermionic degrees of freedom. Assuming the presence of 
just a single species, (spinless fermions) we may describe these with the complex fields 
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and ijjR. Their Hamiltonian contains the standard free part, 



Ho = vfJ dx {ipi{id^ + kF)^L + ^i{id^ - kF)^R) (4.15) 

where vp the Fermi velocity. In Luttinger hquid theory, all (marginal) four- Fermi interac- 
tions are included in the Hamiltonian as well, 

i^int = 27r^f j dx {g2 (V'IV'l) {i^Un) + 54 {Ai^if + (V^^k)') (4-16) 

where g2 and are dimensionless coupUngs. The Luttinger model considered here is left- 
right symmetric, but not Lorentz invariant. It may be solved by bosonization to a single free 
bosonic excitation, whose velocity is renormalized to. 



v = VF^J{l + g^)^-gl (4.17) 

Qualitatively, this effect is similar to the change in velocity we have found for the scalar exci- 
tations in the presence of the charged magnetic brane. There, however, the renormalization 
of velocity affects only one chirality, not both. The renormalization of velocity also affects 
the entropy density dependence on T, as the specific heat coefficient for the above Luttinger 
liquid is given by. 

Comparing this expression with the entropy density for the charged magnetic brane, we see 
that a decreased velocity leads to an increase in specific heat coefficient, an effect that is 
qualitatively consistent with the gravity picture. 
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5 Structure of current and stress tensor correlators 



In this section, we will derive the general structure of the 2-point functions involving the 
U(l) current J and the stress tensor T in the state dual to the charged magnetic brane 
solution. We will show that the correlators may be organized in terms of a Taylor expansion 
in powers of the charge density parameter q of the background solution (defined in (2.17) 
and related formulas). This general structure will be confirmed by explicit calculations in 
subsequent sections. 



5.1 Review of correlators in the purely magnetic solution 

The 2-point correlators of J and T in the state dual to the purely magnetic background 
solution (for which g = 0) were evaluated to leading order in the approximation of small 
momenta The results are given as follows]^ 



(r++(p)r++ 
(r-(p)r-_ 



-p)) 



lg=0 



g=0 



g=0 



kc p+ 
27r p_ 
c 



pI 



487r p_ 

c 
487r 



(5.1) 



All other correlators vanish for g = in this approximation. In fact, since the purely magnetic 
background solution is invariant under charge conjugation, the mixed correlators {JT) in 
the dual state vanish exactly. The overall strength of the correlators is governed by. 



which is the Brown-Henneaux central charge for AdSs with radius 



^-1 



5.2 Extended symmetries for the charged magnetic brane solution 

As q is turned on, the charge conjugation and 2-dimensional Lorentz symmetries of the purely 
magnetic background are spontaneously broken by the non-zero charge density q. Charge 
conjugation reverses the sign of g, while a Lorentz boost transforms g as a (light-like) Lorentz 
vector. The latter motivates the notation, 

q = q+ (5.3) 

'''The properly normalized 2-dimensional current and stress tensor operators, which were denoted by J 
and T in [5], (see equations (4.23) and (6.18) of [5^ for their precise definition) will be denoted respectively 
by J and T respectively here, i.e. without hats. 



25 



where the subscript + stands for a Lorentz index. Extended transformations, which trans- 
form the fields as well as the symmetry breaking parameter, will leave the background 
solution invariant, as always with spontaneous symmetry breaking. Concretely, transform- 
ing all Lorentz vector and tensor fields as well as the parameter g+ under extended Lorentz 
transformations leaves the charged magnetic brane solution invariant. Similarly, the solution 
is invariant under extended charge conjugation. In the dual theory, these extended Lorentz 
and charge conjugation symmetries will be realized through Ward identities for the corre- 
sponding spontaneously broken symmetries. It is this invariance that we will use to organize 
the g_|_-dependence of all correlators for the charged magnetic brane solution. 

Extended charge conjugation symmetry allows the correlators {JT) to be non- vanishing 
functions which are odd in g+ and permits g+-dependent corrections to the correlators {J J) 
and {TT) which are even in g+. Extended Lorentz invariance further restricts the correlators 
as follows. In momentum space, the correlators may be organized through their dependence 
on the combinations q+P- and = p+P- which are invariant under extended Lorentz sym- 
metry, and respectively odd and even under extended charge conjugation. (In particular. 



these are the only combinations that will enter the full reduced field equations (C.5), as will 
be discussed later.) 

5.3 The approximation 

The full analytic calculation of the correlators for general values of momenta is at present out 
of reach, and we will instead construct an approximate solution with the help of overlapping 
expansions. Our aim is to extract the leading behavior of certain correlators in the limit of 
low momenta. It will be helpful to spell out here precisely what this approximation entails. 
In momentum space, the approximation consists in neglecting the following contributions: 

1. purely local terms produced by polynomials in p+ and p_; 

2. non-integer powers of p"^, such as p^'^ and which arise from the near region solution; 

3. positive integer powers of p"^ which arise from the far region solution; 

4. positive integer powers of q+P- which arise from the far region solution for q^ ^ 0. 

Any term which is suppressed by a contribution of the above type relative to a leading term 
will be neglected. As we will see, these rules will allow us to isolate the leading long-distance 
behavior of certain correlators, and will amount to setting other correlators to zero. Within 



this approximation, the correlators of (5.1) in the purely magnetic background solution 
were derived in [6J by explicit calculation using the overlapping expansion methods. The 
corrections to ( |5.1 ) referred to in point 4. are absent here, while the corrections of point 3. 



are immaterial since their presence would produce only purely local terms. 
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5.4 General functional dependence on p± and and analyticity 

Extended Lorentz symmetry restricts each 2-point correlator to be of the form, 

{p+y+{p-y~F{q+p^,p+p^) (5.4) 

where s± are non-negative integers, s+ — s_ is the total spin of the correlator, and F is a 
function of the extended Lorentz invariants q+P- and p"^ = p+p~- Extended charge con- 
jugation symmetry prohibits fractional powers of since these would produce non-trivial 
monodromy under g+ — )■ — g+. The g+ — )• limit will be assumed to be smooth, since it 
physically corresponds to letting the charge density tend to zero, which is a smooth process. 
As a result, F must have a Taylor series expansion in powers of {q+p^)- 



In view of point 2. of the approximations made in section [5^ we are to omit all fractional 
powers of p^, leaving us to retain only integer powers of p^. In view of points 1. and 3. of 
section 5.3, and the fact that s± > 0, we may omit all non-negative integer powers of p"^. 



Putting all together, we find that F must take the form, 

oo 

F(g+p_,J9+J9_) = E E Fr„,n(g+P-)"(P+P-)-" (5.5) 
m=0 n>0 

The sum over m is restricted to even (resp. odd) integers for correlators that are even (resp. 
odd) under extended charge conjugation. Finally, we retain only the leading long distance 



contribution, corresponding to the term in (5.5) with smallest m and largest n. 



5.5 The trace of the stress tensor 



The trace T of the full 4-dimensional stress tensor Tfj,u is governed by the trace anomaly and 
is thus a local function of the metric and gauge field. As a result, any 2-point correlator 



involving T is local and, within the approximations spelled out in section |5.3[ effectively 
vanishes, a result we will denote by T ~ 0. Concretely, T is given ^"'^ 



cv 



K 



(5.6) 



where we have used the form of the conformal boundary metric g^y for the charged magnetic 
brane in the oblique coordinate system (2.5), 



Cv 



dx dx + /3{dx ) + —dx^dx 



(5.7) 



'^Here, and throughout, we set Tii = 722 = Tk, and f3 — —ac/{4:b) with Oj, given in (2.13). By construc- 
tion, the holographic stress tensor is symmetric, so that we have T-\ — T-+- 



27 



Furthermore, it will be proven in section 1J2_ that the 2-point correlator of Tk with any one 
of the operators 7±±, Tj, 
that we also have Tk 



'+-) 



and Tk vanishes within the approximations of section |5.3 
0. Combining this result with T 



so 



and (5.6), we find that 



r+- - fiT- 







(5i 



since its 2-point correlator with any one of the operators J±^ T±±, T+-, and Tk vanishes 



within the approximations of section 5.3 Result (5.8) should come as no surprise: putting 



the metric ds"^ in standard form by changing coordinates to. 



X 



X 



X 



X 



(5x^ 



(5.9) 



condition (5.8) translates to the customary vanishing trace condition 7+- ~ 0. 



The trace condition (5.8 ) leads to relations between the correlators for various components 



of the stress tensor. Some care is needed, however, in keeping these relations consistent with 
the approximations of section 5.3 For example, to order g° the correlators of 71- with 7+_ 

and 7++ vanish, and thus it follows from (5.8) that the correlators of 7h with 7+- and 

7++ vanish at least to order g^, all within the approximations of section 5.3 The correlator 



of 7+- with 71- to order is proportional to f3p^/p+. Contrarily to naive appearances, 
this correlator vanishes within the approximation of section |5.3[ Indeed, since we have 
/3 = —ac/b ~ g^, the correlator exhibits a factor of (g+p-)^, and must be set to zero. 
Analogous arguments hold for the correlators of 7+- with the currents, and lead to their 
vanishing. Collecting the results obtained in this way, we have. 



(r+-(p)r±±(-p)) ~ 
(r+-(p)r+-(-p)) ~ 



(5.10) 



Higher order corrections in will be accompanied by higher powers in p_ and are to be 
neglected as well. Thus (5.10), and thus the operator relation 7+- ~ 0, will be valid to all 



orders in g+ within the approximations of section 5.3 



5.6 The current component J71 



It will be shown by explicit calculation in section 7^ that the following relation holds in all 
2-point correlators within the approximations of the overlapping expansion method. 



J- - Q+T^ 







kac 
CyCE 



(5.11) 



This relation is analogous to (5.8). As a result, the fate of the 2-point correlators involving 
J71, under the approximations of section |5.3| is analogous to the fate of the correlators 
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involving the component of the stress tensor, and we deduce the relations, 



{J-{v)J±{-v)) ~ 

{J-{p)T+-{-p)) ~ (5.12) 
Equivalently, the operator relation J71 ~ will be valid to all orders in ^4., within the 



approximations of section 5.3 



5.7 Conservation of current and stress tensor 

Gauge invariance and translation invariance in guarantee that the current and the 
stress tensor 7^^ must be conserved]^ 

d+J^ + d^J+ -2/3 9_jr_ = 

d+r-+ + d-r++-2(3d-r-+ = o 

d+T— + d^r+- - 2/3 d^T— = (5.13) 

Throughout we set to zero the momenta pi = 0, so that no derivates di occur and Tk does not 
enter. In view of the relations 7+- ~ and ~ derived in the preceding two subsection, 
these conservation relations reduce to chiral conservation equations, 

9_r++ ^ 

d+T— ^ (5.14) 

These equations hold when inserted into any 2-point correlator, and within the approxima- 
tions of section 15. 3[ 



The equations 7+- ~ 0, J7L ~ 0, together with the chiral conservation equations imply 
the general structure of the remaining 2-point correlators. To see how this works, we consider 



first the case q+ = 0, and use these equations to derive the general structure of (5.1). For 
example, by Lorentz invariance, the 2-point correlator of J7+ must be of the form, 

{Mp)M~p))L, - ^/(/) (5.15) 

where / depends only on the Lorentz invariant = p+P-- Following the approximations 
of point 2. of section 5.3, we neglect all fractional powers of p"^, so that / has a Laurent 



^More precisely, the divergence of the current is nonzero in general due to the chiral anomaly, but since 
the anomaly is purely local, it may be set to zero consistently within the present framework. 
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expansion in positive and negative integer powers of p^. Strictly positive powers produce local 
terms, which are to be discarded following point 1. Strictly negative powers are prohibited 
by the chiral conservation equation p_J'+{p) = of (5.14), which implies that / must be 



constant. The value of the constant is derived by explicit calculation [B]. The same argument 

laining ( 

p+T--{p) = within the approximations of section 5.3 



may be adapted to the remaining correlators of (5.1), using the chiral conservation of the 
-ip) ~ 



stress tensor J9_7+ 



5.8 Structure of correlators for the charged magnetic brane 

We will now derive the general structure of the 2-point correlators of and T, constrained 
by extended charge conjugation and extended Lorentz symmetry, in the approximation of 
section 5.3, and under the (mild) assumption that the g+ — )■ limit is smooth and leads to 



the general form (5.5) for all 2-point correlators. 

The relations 7+- ~ and J71 ~ lead to the vanishing of all 2-point functions involving 



these operators, as was already expressed in (5.10) and (5.12). The only remaining non-chiral 
correlator is (7++(p)71-(— p)); its general form is given by (5.4) and (5.5) with s+ = s_ = 0. 
The chiral conservation equations ~ P-T++ ~ readily exclude all terms with n > 1, 

so that this correlator must vanish, 



{T++{p)T-{-p)) ^ Q 



(5.16) 



The remaining correlators are all chiral. 

For negative chirality, only a single 2-point correlator, namely {T--{p)T- 
in view of JTI ~ and (5.12). Its general structure is given by (5.4) and (5.5) with s+ = and 



remains 



4. Chiral conservation p+T- ~ forces us to restrict to n = 1. On the other hand, all 



terms with m > are to be omitted by point 4. of the approximations of section 5.3 Hence 
only the m = 0, n = 1 term survives, which is independent of g+ and thus must coincide 
with the result from the purely magnetic case. 



{r-ip)r-i-P)) 



p. 



487r p+ 



(5.17) 



For positive chirality, we find the following 2-point correlators, 

{J4p)J4-p)) 
{j+{p)r++{-p)) ~ 71- 



kc 

271 p_ 
Q+P+ 



P- 



487r j9_ 



P- 



30 



where 71 and 72 are two coefficients which are not determined by general arguments, and 
which will have to be obtained by explicit calculation. The terms proportional to c on the 
first and last lines correspond to the g+ = contribution from (5.1). 

To prove the form of the correlator {J+J^) in (5.18), we set s+ = 2 and s_ = in (5.4), 
and use chiral conservation p_J'^{p) to restrict the expansion of F in (5.5) to n < 2 
and n < m + 1. The n = 1 contribution is local for m > 1, while the n = 2 contribution is 
always localj^ Hence only the n = 1, m = contribution remains, which is independent of 



g+ and coincides with the ffist line of (5.1) which yields the first line of (5.4). 

The form of the correlator (J7+7++) in (5.18) is proven in analogous fashion. Its general 
form is given by (5.4) and ( 5.5[ ) with s+ = 3, s_ = 0. Chiral conservation of either J7+ or 7++ 
forces n < 3 and n < m + 1. Extended charge conjugation requires m to be odd, so that we 
must actually have m > 1 and n > 2. The contribution n = 3 is local, leaving only n = 2 
with m = 1, with all m > 3 contributions to be omitted in view of point 4. of section 5.3 



This is precisely the form given in the second line of (5.18). 

The form of the correlator (7++7++) in (5.18) is given by (5.4) and (5.5) with s+ = 4 and 
s_ = 0. Chiral conservation of 7++ forces n < 4 and n < m+1. Extended charge conjugation 
requires m to be even, which restricts to m = 0, 2. All contributions with m — n > 1 are 
to be omitted in view of point 4. of section 5^, while the n = 4 contributions, as well as 
the n = m = 2 contributions are local. This leaves only the contribution m = 0, n = 1 

2, = 3 
its 



corresponding to the first term of the correlator in (5.18), and the contribution m 



corresponding to the second term in (5.18 



). Since the first term is independent of q 



normalization is provided by the result of (5.1). 



Explicit evaluation in section [8] will con&m the general form of the 2-point correlators 
derived above, and will determine the values of the coefficients as follows. 



71 



CvCe 



4:nG3q+ 



72 



„2 2 



4:7TGskb{q^ 



(5.19) 



The coefficients are independent of g+ since CE/q+ is. We will also show in section Is] that 



the 2-point correlators in (5.18) remain the dominant IR contribution throughout k > 1/2. 



5.9 Structure for general value of B > Be and twisting 

The value of the normalized magnetic field B, which was expressed in terms of the parameters 



of the charged background solution in (2.18), is found to enter the general structure of the 



2-point correlator only through 71 and 72. In particular, for the purely magnetic solution 
we have B = 00 and g+ = 0, so that 71 = 72 = 0, but for finite B we have 71, 72 7^ 0. Thus, 
the two Virasoro algebras and the chiral U(l)-current algebra, which were identified for the 



^Note that a contribution in (p_) " with n > 1 is local, even though it is not polynomial in p±. 
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purely magnetic case, persist for all B > Be- As the charge density is turned on, however, 
the current and stress tensor operators become twisted. In the positive chirality sector, this 
may be seen by diagonalizing the system of 2-point functions of J'+ and 7++, by setting. 



(0) 



/++ — fJ^+J- 



(0) 



CyCE 

kb 



27r 
kc 



(5.20) 



If and 7^? obey the 2-point functions of the purely magnetic solution of (|5.1l), then J7V 



and 7++ defined by (5.20) will obey (5.18). Note that this result crucially depends upon the 



precise value taken by 72, and its relation to 71, given by, 

27r o 



72 



kc 



7i 



(5.21: 



Note that the twisting of (|5.20[) is possible because both j7+°'' and 7++ are chirally conserved. 



-(0) 
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6 Linearized metric and gauge field fiuctuations 



In this section, we begin the process of calculating the two-point correlators of the U (1) cur- 
rent and the stress tensor T. To do so, the full Einstein-Maxwell-Chern-Simons equations 
are linearized around the charged magnetic background solution. Identifying the source and 
vev components in the Fefferman- Graham expansion of these solutions, we derive the vevs 
of the current and stress tensor using (2.4), and from there extract (the non-local parts of) 
all two-point functions. Since the background solution is invariant under translations in x'^, 
we solve for plane waves with fixed momenta for fi = +,—, 1, 2. For simplicity, we restrict 
to the most interesting case where pi = P2 = 0, leaving only the momenta p±. 



6.1 Parametrization of linear fluctuations 

A general Ansatz for the plane wave fluctuations of the metric and gauge field is given by. 



A = Ac + ttMe'^'^dx^^ 
ds^ = dsl + hMNe'P^dx^dx^ 



(6.1) 



We continue to use the notation px = p^x^ + p^x . Here, Ac and ds^ are the gauge 



potential and metric of the charged magnetic brane solution. In the gauge adopted in (2.5) 
referred to as LMN gauge, the charged background solution takes the form. 



dsl 



dr"^ 



+ Mo{r){dx+f + 2Lo{r)dx+dx- + e^^''^'''^dx'dx' 



A 



c 



hx^dx^ + Ao(r)c?x^ 



(6.2) 



where the functions Lq and Vq are those of the purely magnetic background solution, obeying 
equations El, E4, and fV of (2.6) with = P = 0, while Aq and Mq are given in (2.11). 



Translation invariance in x^ of the solution (6.2) justifies our consideration of plane wave 



fluctuations in (6.1). The functions aM and qmn depend on the holographic radius, but are 



independent of x^. Invariance under rotations in the x^'^ plane of the solution (6.2) leads 



us to require rotation invariance of the fiuctuations of (6.1), which implies the following 
conditions, 



hi. 







hij 



Siih 



ijlLK 



1,2 



(6.3) 



Even with the benefit of translation and rotation symmetry, the corresponding reduced field 
equations remain quite involved, and will not be presented here. 

To obtain workable reduced equations we choose a gauge. It turns out that a convenient 
gauge choice is not global, but is rather obtained by patching together the light-cone gauge 
for the near region and the LMN gauge for the far region. In the overlap region, a gauge 
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transformation is required to patch the two gauge choices together. In the extreme UV hmit 
of the far region, we will make a further gauge transformation to Fefferman- Graham gauge, 
which is required to properly extract the current and stress tensor data. 



6.1.1 Near-region: light-cone gauge 

In the near region, defined by r ^ 1, we will adopt light-cone gauge, and denote the corre- 
sponding fluctuation fields with tildes, namely cim and Hmn- These fields satisfy the usual 
light-cone gauge conditions, 

a_ = hr- = = h__ = (6.4) 
in addition to the consequences of invariance under translations in x^, and rotations in x^, 



as in (6.3). The full gauge field and metric in hght-cone gauge become, 

A = Ac+ {hrdr + a+(ix+) e*^^ 
ds^ = dsl; + {h„dr'^ + 2hr+drdx^ + h++{dx^f + hKdx^dx'^ e^P"" (6.5) 

where the remaining fields hrr, hr+, and hx depend only on r. 

6.1.2 Far region: LMN gauge 

In the far region all momentum-dependence in the reduced field equations is to be neglected. 
This specification always requires p+p- ^ r. As we have discussed for the case of scalar 
fluctuations in section 3.2, however, depending on the value of k, additional restrictions on r 
may be needed in order to ensure that the far region will extend non-trivially into the near 
region. We will exhibit these conditions as we proceed forward. We stress, however, that 
even if the reduced fleld equations do not involve momenta in the far region, their boundary 
conditions will be given by p-dependent coefficients. 

When all momentum dependence in the reduced equations is neglected, the fluctuations 
will be analyzed in LMN gauge, where they correspond to x'^-independent fluctuations of the 



functions L, M, N,V, A± of (2.5). LMN gauge will be convenient here, because this linear 
fluctuation problem was considered and solved already in Section 5 and Appendix B of [H] . 
The fluctuation flelds in LMN gauge will be denoted by the original letters om and Hmn, 
and are subject to the following gauge conditions, 

2 . , Mo, 
73^+- + JA 



hr+ = = hrr = — 7-3/?-+- + ^4:^ (6.6) 



in addition to the consequences of invariance under translations in x'^, and rotation in x^, 



X 



2 



as in (6.3). The relation for hrr is imposed to guarantee that the form of the dr term is 
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preserved by the fluctuations. Thus, the gauge fleld and metric in LMN gauge become, 

A = Ac+ (a+dx"^ + a_(ix~) e'^'' (6.7) 
ds"^ = dsc + {hrrdr"^ + h+^^dx'^f + 2h+_dx'^dx' + h__{dx~f + hxdx'dx'^e'P'' 
The remaining functions a±, hrr, h±±, and hx depend on r, but are independent of x^. 



6.2 Reduced field equations in the near region 

In the near region, where r ^ 1, we adopt hght-cone gauge. The reduced equations for 
fluctuations in hght-cone gauge around the full charged magnetic background solution are 
derived in Appendix B. For r <^ 1, it is appropriate to use the approximation! 



Voir) 
Loir) 



-- 
2br 



Aoir) 
Mo(r) 



qr 

IT 



2q''r 



ki2k 



(6.8) 



so that the general light-cone equations (C.l), (C.2), (C.3), (C.4), and especially (C.5) 



greatly simplify. The reduced equations may be presented most symmetrical in terms of the 
functions hx and c(r) = 4ra(r), and we have. 







r^h';, + 2rh'^-hx-J^^h 



K 



^hx -\ 77, Cr 



72ki2k 



= r'c; + 2rc, - kik - 1)5, 



-Cr 



9 9 9i"— 9 

q^p^r iqp^r 
-Cr - 



k-l 



72ki2k 



63 



-h 



K 



(6.9) 



Note that the equations only depend on the combinations p^ = p+P- and For later use, 
we also record the relations giving the remaining flelds in the near region. 



Qj.2 

4ibr ~ , 
h 



r hj^ + 2rh^ — -hx + T^qp~r a, 



P- 

Rescaling r by p^. 



K 



pt 
24:ibr'^ 

P- 



63 



iirttr)' + kar 



p^x 



hxir) = hxix) 



Cr(r) = Crix) 



(6.10) 



(6.11) 



^In keeping with the discussion of section 



2.7 



we restrict attention to the critical background solution for 
which B — Be, and ag = 0; so that no linear term in Mq is present as r — >■ 0. The solutions for B > B^ may 
then be reconstructed by applying an SL(2, R) transformation to the critical solution. 
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reveals that the equations in terms of the rescaled variable x, 



Ax 







X c" + 2xc^ — k{k — l)Cr 



Ax 



(6.12) 



intrinsically only depend on the combination which was already introduced in (3.23) for 
the scalar fluctuation problem. Here, we have used the notation = 2k{2k — l)/3. 

On the one hand, for > 1, we automatically have ^ — )■ in the limit of small momenta 



p±, no matter how this limit is being taken. Also, the .^-dependent terms in (6.12) are 



well-behaved as x 0. On the other hand, for 1/2 < k < 1, the combination ^ exhibits 
a directional singularity in the limit of small momenta p±, so that ^ may tend to zero in 
one limit, to oo in another, or to any finite value. Also, for 1/2 < A; < 1, the ^-dependent 



terms in (6.12) diverge as x — )■ 0. Finally, it is unlikely that equations (6.12) can be solved 



for finite ^, beyond the type of perturbative results that were derived for the scalar problem 
in Appendix A. To avoid these three complications, we will henceforth restrict to the case 
where A; > 1, and proceed to solve it in the limit of low momenta, which implies ^ ^ 1. 



6.3 Solutions in the near region 



Assuming k > 1 and ^ ^ 1, equations (6.15) decouple from one another, 








Ar 



r Cj. + 2rc,. — kik — l)cr — —c. 

Ar 



(6.13) 



and may be solved in terms of modified Bessel functions. 



hK{r) 



CAT] 



Co 



'2cr+l 



'2fe-l 




(6.14) 



where vq and cq are independent of r, and a was given in (2.16 ). As was the case for the scalar 



fluctuations, the above linear combinations of Bessel functions are required by regularity of 
h and a as r — )■ 0, and are proportional to K2a+i and K2k-i respectively. 
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6.4 Solutions in the overlap region (light-cone gauge) 



The overlap region, defined by ^ r ^ 1, is contained in the near region so that the 



approximations of (6.8), as well as the reduced equations of (6.13) are valid. It is also 
contained in the far region <^ r, where all momentum dependence in the reduced field 



equations is to be neglected, so that equations (6.13) actually reduce to 



Its solutions are given by. 



= r^h';^ + 2rh'^-hK 



(6.15) 



(6.16) 



The coefficients v± and c± are independent of r, and may be determined in terms of Vq and 
Co respectively by matching the solutions of (6.16) with those of (6.14) in the <^r limit. 



Here, we will only need the ratios of these coefficients 
Ti-2a] 



r(2 + 2o-) V2 



C_ 



r(2-2A;) fp\^^-'^ 
.2. 



V{2k) 



(6.17) 



Clearly, to leading order in small momenta, the coefficients f _ and c_ are suppressed com- 
pared to f+ and c+, and the corresponding terms may effectively be dropped. In terms of 
the coefficients f±, c±, the remaining light-cone gauge functions are given as follows. 



hK{r) 

g^2 



P- 



av,r'^ — (a + l)v^r 



-cr-l 



Qib(2k - 1) . 

a+ = — ^ -c+r'^ 

P- 



(6.18) 



The cancellation of is brought abou t by the fact that its various terms in (6.10) cancel 
for the overlap region solution of Hk of (6.16) and/or are of higher order in p'^/r. 



6.5 Change to LMN gauge in overlap region 

We have solved the linearized equations around the charged magnetic brane in the near 
region in light-cone gauge, and need to match the result with the solution in the far region 
in LMN gauge. The matching takes place in the overlap region ^ r ^ 1 where both 
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approximations are simultaneously valid, and where a simplified solution holds. The general 
form of the solution in light-cone gauge will be denoted by, 



Ac + dM{r) e'^'rfx^ 
HMNdx^'dx'^ = dsl + hMN{r)e'P''dx^dx^ 

while its general form in LMN gauge will be denoted by, 

- Ac + aMir)e'P''dx^ 

dsl + hMN{r)e'P''dx^dx^ 



AMdx 
Hiuwdx^^ dx^ 



(6.19) 



(6.20) 



Since the background solution is the same in both gauges, the gauge transformation be- 
tween the linear fluctuations may be carried out at the linearized level. This involves a 
diffeomorphism vector field U^^ and a gauge transformation field G which take the form. 



U^\r,x) 
6(r, x) 



e{r) e'P" 



(6.21) 



To convert the metric and gauge fluctuations Hmn and Am of the light-cone gauge to Hmn 
and Am of LMN gauge in the overlap region, we solve the equations. 



Hmn + ^ mUn + nUm 



Am + U'^Fkm + dMQ 



(6.22) 



Hmn 
Am 

for ^ 9, as well as Hmn and Am- The covariant derivatives V m, and the field strength 
Fkm are with respect to the background solution, and so is the metric used to raise and 



lower indices. Substituting the plane wave forms of (6.19) and (6.20) 



/iMTve'P^ + VmUn + VnUm 
fiMe^P" + U^'Fkm + BmQ 



(6.23) 



Imposing rotation invariance, as spelled out in (6.3), on hMN, clm, h MN, and qm gives 
^ — "'^ — 0. Next, we impose the conditions for light-cone gauge of (6.4) on Hmn and 



u 



u 



qm, and for LMN gauge of (6.6) on Hmn, and qm- Finally, we assign the expressions for 



the light-cone gauge solution in the overlap region of (6.16) and (6.18) to t he re m ainin g 
components of the fields Hmn and qm- The resulting equations are derived in (D.2), (D.3), 



and (D.4) of Appendix D 



6.6 Solutions in the overlap region (LMN gauge) 

The corresponding solutions in the overlap region for the metric and gauge fluctuations in 
LMN gauge are constructed in Appendix [D] as well, and are of the following form, 

h++ir) = s++r + t++ 
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with the coefficients given by, 

s = Aihp-UQ 

s+„ = 

(T_ = IP-Oq 



h (r) 

h+-{r) 
a+(r) 
a_(r) 



s r + t_ 



r 1 r 



(6.24) 



+ 



r_ 



246 ' 

246 
= 2hul + 

(2A; - 1 
Ak 



+ 



P4 



246p_ 

,3 



-S 



pi 



2Ahp+ 



1 

246 

6z6c. 

I 

P- 



2 I 2 

p.s + p_s_ 



Abp^k Abp^ 



(6.25) 



Since c^/c^ ~ , we may effectively set c_ = within the approximations of section 5.3 



The metric fluctuations of the field are given by, 

hK{r) = v+r" + v-r-^-" (6.26) 
and decouple from all other fluctuations of both the metric and the gauge field. Since 



V-/v^ ~ p^"""*"^, we may effectively set f _ = within the approximations of section 5.3 
6.7 Reduced field equations in the far region 

In the far region, defined by p^ <^ r, all momentum dependence in the field equations is to be 
omitted. The corresponding reduced field equations for the metric and gauge field fluctua- 
tions, in LMN gauge, coincide with the perturbative expansion of the fields L, M, N, V, E, P, 
obtained already in equations (5.35), (5.36) and (5.37) of [H]. Making the correspondence 
with the present notations, the expansion is as follows. 



L(r, 


x) 


= Lo(r)+£Li(r)e*J'" 


Li{r) = h^_{r) 


M(r, 


x) 


= Mo{r) + eMi{r) e'P"" 


Mi(r) = h++{r) 


N{r, 


x) 




Ni{r) = h^^{r) 


E{r, 


x) 


= Eoir)+eE,{r)e'P'' 


Eiir) = a^ij-)' 


P{r, 


x) 




Pi(r) = — a_(r)' 


V{r, 


x) 


= Voir) + eV.ir) e'P'' 


2e2^°("Vi(r) = hxir) 


fir, 


x) 


= /o(r)+s/i(r)e^P^ 





(6.27) 
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where e is a formal expansion parameter. In [TT] , we had also introduced the gauge potentials 
Ai and Ci, which are related to the fields Ei and Pi by Ei = A[ and Pi = C[; these potentials 
are related to the present fields by 



-Ci{r) 



a+{r)-Al 
a_(r) - A° 



(6.28) 



where A^ are constants. Here, / was defined in (2.7), and we have. 



/o 
/i 



^0 



2LoLi - MqA^i (6.29) 
Maxwell's equations for the perturbation functions are given by. 
Ml = [e^^''{EoNi + LqPi))' + 2A;6Pi 

M2 = (e^^o(2LoPo'^4 + ^o^i + ^o^i + MqPi))' - 2kbEi (6.30) 

while Einstein's equations are given by. 

El = L'l + 2V^L[ + 2VlL'o + 4Li (fJ' + (V^f) + 4Lo (K' + 2V^Vl) - AE^Pi 

E2 = M'l + 2V^M[ + 2V;'Mo + 4Mi (v^;' + [V^f) + 4Mo {v' + 2V^oVi') + SE^Ei 

E3 = N'l + 2V;,N[ + ANi (vj' + (^J)') 

fV = + 4Fo7( + 4V;7o + 2K7o + 2^0"/! + 4(K,')'/i + S^IK/o 

E4 = (6V^o" + 12(K,')') /i + (6^/' + 24FoV/) h + 6^//^ + 6Vo7{ 



-?,2h'^e-^^''Vi + SLoEqPi + AN^El 
We also have first integral equations, of which we will need the following. 



(6.31) 



e^^»(EoiVi + LoPi) 
e2^°(2Lo^o^i + L^Ei + Po^i + MqPi) 
(A^iMo - MoA^Oe^^" - ^khAoCi 



-2kbCi 

2kbAi 

2Ao 



(6.32) 



The solution to these equations in the far region needs to be matched to the solution of 
(6.24), (6.25), and (6.26) that we have derived earlier for the overlap region (recall that the 
overlap region is contained in the far region, so that both solutions must match in the overlap 
region). Thus, in the overlap region, we must have. 



Li(r) 



-r + t. 
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s r + t_ 



Mi(r) : 
iVi(r) : 
Ai(r) : 
-Ci(r) : 

where the coefficients s±±,t±±,S4._ 
The coefficients v± are decoupled from all other fields. At low momenta, the coefficient V- 
may effectively be set to 0. 



+ f _r" 



(6.33) 



,t+_,cr-|-, and T± are related by the equations of (6.25) 



6.8 Solutions in the far region 



The system of linear equations of (6.30) and (6.31) for the far region turns out to be solvable 



by quadratures in the following sense. Their general solution may be obtained by (successive) 
quadratures in terms of the functions Lq and Vq which characterize the purely magnetic 
solution. In fact, Lq itself may be obtained in terms of Vq by quadratures. Thus, although Vq 
is not (as of yet) known analytically, all other fluctuation functions are determined explicitly 
in terms of Vq. The system may be integrated iteratively in the following sequence, 

E3 — ^ Ml — y (fV,E4) — ^ M2 — ^ E2 

Ni ^ Pi ^ ifuVi) ^ El ^ Ah (6.34) 

The top line gives the order in which the equations should be solved, and the bottom line 
gives the corresponding functions obtained. Equations fV and E4 are coupled and must 
be solved together; their solution is obtained by solving a single auxiliary third order linear 
differential to which two solutions are a priori known for symmetry reasons. The full solution 
in the far region was obtained in |llj . 

In preparation for the calculation of correlators, we will review here the solution for the 
functions A^^i and Pi, which do not require the more complicated solution for Vi and /i. 



These solutions were obtained already in [TT]. Inspection of equation E3 for A^i in (6.31) and 
equation El of (2.6) shows that one solution for A^^i is proportional to Lq, while the other 



(linearly independent) solution is proportional to the function Lq, deflned by, 

dr' 



Loir) 



Loir) f 

J CO 



Lo(r')2e2^«M 



(6.35) 



The function Lq involves data of the purely magnetic background solution only, and has the 
following asymptotics. 
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(6.36) 
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Thus, the general solution for A^i is a linear combination of Lq and Lq. Since we have 
Lo(r) = 2br for small r, we see that the asymptotic behavior of Ni{r) is indeed of the form 



given already in (6.33) for the overlap region. Identifying the coefficients gives, 



iVi(r) 



2b 



2bLl{r)t. 



(6.37) 



Since this solution does not involve the charge density parameter g, it coincides with the 
corresponding fluctuation solution for the purely magnetic case of [6]. 

Next, the function Pi, and thus A_, may be obtained by solving equation Ml. It is 
actually more convenient to solve instead the first integral equation on the third line of 
(6.32), since this equation gives directly Ci = — A_ in terms of data that we have 



already obtained. The result is as follows. 



A. 



(J- + 



Aq 

AkbAQ{r) 



+ 



Mr) 
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s 



Mo 
AkLnA 



(6.38) 



0^0 



Using AQ{r) ~ qr^/k as r — )■ 0, we match the asymptotics of this expression with the one in 



the overlap region given in (|6.24|). This accounts for the presence of the additive constant 

p'^^ is suppressed 



(T_ in (6.38), and solves for Aq = ^hqr^. Since the coefficient r_ ~ c_ 



within our approximation, we will set 



r_ = Ao = 



(6.39) 



in (6.38), and throughout the remainder of this paper. 
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7 Current and stress tensor correlators: general set-up 



In this section, we will begin the calculation of the actual 2-point correlators of the current J 
and the stress tensor T in the state dual to the charged magnetic solution. A number of 
general implications will be derived from the structure of the sources and the solutions in 
the far and overlap regions, which will form the basis for the computation of all correlators 
in later sections. 



7.1 Identification of sources and vevs 

Sources and vevs may be identified from the r — > oo boundary asymptotics of the metric 
and gauge fluctuations, following (2.2) and (2.3), 
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iVi(r) 
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rsj 


Af + ^Af + ... 


-C,[t) 


rsj 


4r 
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1,2 



(7.i: 



The factors of 4 in the normalization of r have been included to ensure that the conformal 
boundary metric has the proper customary normalization, given by, 

a 

^ ~h 

which requires the following relation between r and the Fefferman- Graham coordinate p. 



dx^dx~ + f}{dx'^Y + ^dx'dx' 



(7.2) 



P 



. h 3 + 61nr 

4(r - ro) + , + 



(7.3) 



where tq and l\ are constants which will not be need ed h ere. The current and stress tensor 
vevs are given in terms of g'^^} and by equation (2.4). As in |6j, we will prefer to deal 
with normalized 2-dimensional current and stress tensor, obtained by integrating over the x^'^ 
space with coordinate volume V2, and effective 3-dimensional Newton constant G-^ = G5/V2. 
The expressions for the vevs of these normalized current and stress tensor are given by, 

8ttGsJ^{x) = cvA^^\x) + local (7.4) 
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We close this subsection by making a key observation. The sources to the metric and 
gauge Geld fluctuations, gj^J and A^^^ are independent of the background charge density g+. 
Indeed, the sources are external fields whose strength may be dialed from the outside without 
reference to the value of g+. This important observation will be used throughout. 



7.2 Fluctuations in Vi and correlators with Tk 

In this section, we will derive the 2-point correlators involving the stress tensor component 
Tk = Til = T22 conjugate to the metric variations in Vi. Using the general result that the 



trace of the full stress tensor T, defined in (5.6) is governed by the trace anomaly, and has 



local correlators with all operators, we will derive the result of (5.8) for the two-dimensional 
trace part 7+_ — /571_ ~ 0. 

The solutions in which Tk is sourced are governed by the fiuctuation equations for /i, Vi. 
These equations were solved in Appendix B of [11], and the general solution takes the form. 



Vi = (tVl + QV,^ + CnVr + VI (7.5) 

where the scripts t, d, n refer to the three homogeneous solutions for which Ni = Pi = 0, 
while p stands for the inhomogeneous solution arising from Ni,Pi 7^ 0. Without loss of 
generality, the latter may be normalized so that Vf vanishes at r = 00 and r — )■ 0. The 
mode V{^{r) behaves as r~^~'^ for r — )■ 0, matches onto the coefficient f_ in the overlap 



region, and is thus to be omitted within our approximation of section pl3\ so that Cn = 0. 

In summary, the stress tensor Tk is sourced only by the dilation mode V{^{r), so that we 
set Ct = Cn = Ni = Pi = 0. The remaining fiuctuations for the dilation mode are given by, 

Lf{r) = rLo(r)' - Lo(r) 

M{^(r) = rMo(r)' - Mo(r) 

Efir) = rEoir)' -Eo{r)/2 

V^^r) = rVoir)' (7.6) 

This is an exact solution to all orders in g+. Since Ni{r) = 0, and Mf{r) ~ r^'^ as r 0, 
the metric matching conditions in the overlap region imply that we must have, 

s±± = t±± = (7.7) 

Since Lf{r) ~ r^"'"'^, Ef{r) ~ r^~^, and Vi{r) ~ r"" as r — > 0, we may always match these 
solutions of the far region in the overlap region by adjusting Uq and c+ in (6.25), and v+ 



in (6.26) without creating any relations between the sources and the expectation values. 
Finally, since Vi{oo) = 1, the coefficient (d is given by the source gl^^ for i = 1,2. As a 
result, the expectation values of the operators 71 _, 7+_, 7++, J7+, JTL, and Tk are all local 
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functions of g^^\ The corresponding correlators vanish, and by locahty of the total trace T, 
so do the correlators with the combination + PT~^~^, and we have, 
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These results form the basis for the derivation of the relations (5.10). 



(7. 



7.3 Expansion of A^i in powers of 



The solution for the field Ni{r) was constructed in (6.37) 

Loir 



iVi(r) 



2b 



-2hLl{r)t. 



(7.9) 



The functions LQ{r) and i^o(r) are data of the purely magnetic background solution, and 

thus independent of g+. Therefore, all dependence is contained in the coefficients s and 

t , which themselves admit a Taylor expansion in g+. In fact, a Taylor expansion holds for 

all components of the tensors Sfj,^ and t^^, and we haver 



M I +[1] 



t 



[2] 



(7.10) 



It is important to notice that, even though A^^i is a metric fluctuation, its expansion in powers 
of g+ may contain odd powers through the gauge fleld sources. 



The relations between s±± and t±±, given on the flrst two lines of (6.25) may similarly 
be expanded in powers of g+. Since the coefficients of these relations are independent of g+, 
the relations are found to hold order by order. 
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p+p- 
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P+P- 
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246p_ 

pt 



'^ 246p+'++ 

Next, we use the key observation of the end of section 



(7.11) 



7.1 



the sources (7^°) and A^f^ are 
independent of g+. In particular, the source of the fleld Ni must be independent of g+. But 



^The order of expansion in powers of q+ will be denoted by a square bracket superscript, while the 
Fefferman- Graham order of expansion in powers of will be denoted by superscript parentheses. 
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this source may be identified from tlie r — > oo asymptotics of the solution by (7.1 ) and, using 
Lo{t) ~ 2br and Lo(r) ~ —1/(26) is found to be given by, 4:g^^l = s As a result we have, 

[0] 1 (0) 

s'!'L = n > 1 (7.12) 



and thus (7.11) simplifies as follows. 



P+P- >] 
246 

P- In 



2Abp+ 



n > 1 
n > 1 



(7.13) 



which implies p^t^L = t++. The r — t- oo asymptotics allows us to extract the expectation 
value of the normalized stress tensor component T and we find, 

167rG3T__ = = 26t__ 

Note that all (even) orders in g+ are allowed to contribute to t 



(7.14) 



7.4 Expansion of in powers of q+ 



The solution for A_ obtained in (6.38), with the overlap region matching conditions of (6.39), 
is given by, 

and also admits a Taylor expansion in powers of g+. Let us now investigate how this expan- 
sion affects the source and vev for this field. To derive the r — >■ cx3 asymptotics of A_, we 
make use of the fact that Lq/Lq = (9(r~^), as well as the following auxiliary asymptotics, 

Mr) ~ ^S^fl-^Uo(r-) 



kb \ Cyr ^ 

+ C(r-2) (7.16) 



Lo(r) 26 

As a result, we find the following equations for the source term A^!^\ and the vev term A^^\ 

/l{0) , CvCiJ , 

A_ = cr_ + — — ;-s + 1 



12k SCyCE 
,(2) _ I «cfc6 

6 2cIce' 



^-^ = -e^- + 73^t- (7-17) 
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Using the fact that the source A_' is independent of g+, and the ear her result (7.12), we 
obtain the foUowing relations between the expansion coefficients in powers of g+, 



A8k ^ ScyCE 

= ^N + ^f^^t"-!] n>2 (7.18) 

- SCyCE - ^ ^ 

where we have used the fact that ce and ac/cE are linear in The expectation value of 
the normalized current component J7L is given by, 

8nG,J. = cyA^'^ = ^t__ (7.19) 

2cyCE 



where we have omitted the local term proportional to s in deriving this result from (7.17). 



7.5 Proportionality of J71 and T 

Comparing the formulas for the expectation values for J7L in (7.19) and for 71_ in (7.14) 
we readily obtain a simple identity between these expectation values. 



Q- 



2cyCE 



(7.20) 



In deriving this result, we have omitted the fractional power dependence on by setting 
c_ = w_ = 0. Since the result is obtained by overlapping expansion methods, we are of course 
also omitting the integer power corrections in and q+p^- In sum, relation (7.20) holds 



within the approximations of section 5.3 Since (7.20) holds for any sources, it is equivalent 
to the operator identity, ~ Q+Tl- announced in (5.11). This equation is invariant under 
extended Lorentz symmetry, since (5+ is linear in g+. 
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8 Calculation of the correlators of Jj^ and 7++ 



Using the general set-up of the preceding section, all correlators involving and 7++ will 
be evaluated in this section. The calculations will be organized along the lines laid out in 
section [Sl and within the confines of the approximations discussed in section|5.3l Correlators 

— 

involving J7+ will be reached by turning on the single source A>_ , while leaving the remaining 



sources and turned off. Similarly, correlators involving 7If+ will be reached by 
turning on only the source g^^_. The chiral correlators (j7+j7+), (J7+7++), and (7++7++) will 
be evaluated in this manner to all orders in g+. We will also take the opportunity to check 
that correlators of Jj^ and 7++ involving mixed chiralities, or the trace of the stress tensor, 
vanish as announced in (5.10) and (5.12). 



,(0) 



8.1 Sourcing J 

We begin by sourcing only the current by turning on y4?\ This case is the simplest, 
and we will discuss the corresponding calculation in detail. The general solution for the 



reduced Maxwell fields a±(r) in the overlap region was derived in (6.25). Using (6.39), it 
may alternatively be expressed as follows, 

a+(r) = (T+ + — Ao(r) 



Here, we have used the r — asymptotics of the function v4o(r) in (2.16) to recast (6.24) in 



terms of in the overlap region. The corresponding Maxwell field strengths are given by. 



Pi{r) 



(1+ 







(8.21 



We will now show that the form of the perturbations exhibited in (8.1) and (8.2), although 



initially derived for, and valid in, the overlap region only, actually provides an exact solution 
valid throughout the far region, under certain conditions. To establish this, we point out 
that the vanishing of Pi, together with the proportionality of Ei{r) to Eo{r), indicates that 
these perturbations may be viewed as the result of transforming by. 



(1+ 



kr. 



(8.3) 



Since r+ should be viewed as a linear perturbation, this shift is infinitesimal. It may be 
interpreted as a rescaling of g+ and may be achieved by an SL{2, R) transformation of 
(2.20) with Ao = 1 + kT-^/q+ and A = 0. Here, r+ actually depends on momenta and on 



source functions, but from the point of view of the differential equations in r for the far 
region, this extra dependence is inconsequential. 



48 



8.1.1 Solution in the far region 

The coefficients a± and r+ may be eliminated in terms of the sources 



field A±. Using the large r asymptotics result for v4o(r) given in (2.13), and the relation 
p+cr„ = p-cr+ derived from (6.25), we find. 



v4i°^ to the Maxwell 
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a_(r) 
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CyCE 



Since the solution (8.4) is effectively generated by the transformation of (8.3), we readily 



determine the solutions in the far region for the remaining functions, and we find. 



Mi(r) 



Li = iVi = Pi = = A 

CyCE 



1.5) 



The first two terms in Mi provide homogeneous solutions to the equation E2 for Mi(r), 
while the third term provides an inhomogeneous solution sourced by the perturbation Ei{r). 



In partial summary, the fields of (8.4) and (8.5) solve the far region equations, and their 



Maxwell fields properly match onto those of the overlap solution. To make them into full 
hedged solutions, it remains to ensure that also the metric functions match with those of the 
solution in the overlap region, and satisfy the absence of metric sources as r — t- oo. 



8.1.2 Matching the metric solutions in the overlap region 

Matching with the overlap region solution determines the coefficients, 
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and insisting on vanishing metric sources as r — )■ oo gives the following relations. 
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■'It will be instructive to temporarily keep both sources A'^'' even though we will ultimately set A';^' — 
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With the help of these equations, we evaluate r^'^ using the overlap region solution of (7.11 ) 
and (7.12) for g^l_ = 0, and we find, 
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2kbac p+P- 
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Clearly, t^^_^ will provide only local terms in the expectation value and may be omitted 
within our approximation. The term sourced by in tt^L similarly produces only local 
terms and may be omitted. Thus, if we set = 0, then all contributions to t±± are local 
and may be omitted in the ma tchin g pro cess in the overlap region. In conclusion, when 
A^^ = 0, the expressions of ( 
match with the solution in t 



(8.4) and (8.5) provide an exact solution to all orders in and 
:he overlap region, up to local terms. 

To close, it is interesting to see what happens when ^ 0. There is now a non-local 
contribution to A^^i, forced upon us by the matching of the solution in the far region with 
the solution in the overlap region. But this contribution, in turn, will sou rce Pi as well 



as Vi, fi, Li, El and Mi. Therefore, when ^ 0, the expressions of (|8.4[) and 



incompatible with the overlap solution, even within the approximations of section 5.3 



8.5) are 



8.1.3 Calculation of correlators with Jj 

under the assumption that A^^^ 



Putting together the result of section 
the following solution. 



8.1.2 



0, we obtain 
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(8.9) 



CyCE P- 

with the remaining fields Li = A'^i = Vi = 0, up to local terms. This solution is valid to all 
orders in g+, within the approximation of section 5.3 Using the r — )■ cxd asymptotics of Ao 
given in (2.13), and the following r — )■ oo asymptotics for Mq, 



Mo(r) + ^Lo(r) 
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kbr 



(8.10) 
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we derive the asymptotics of the fields as r — )• oo, and we find, 
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p_ r 

The expectation values of the current JT"" with all components of and T may be read off 



from ( |8.11[ ) and Li = A^i 



to vanish, J. 

implies that the corresponding correlators also vanish 



Vi = 0. First of all, a number of expectation values are found 
= Tk = 0, within the approximations of section [5^ This result 



[J-ip)J^i-p)] 

[r--{p)J-{-p)] 
[%-{p) J- {-P)] 
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Converting the upper index to lower indices, 

J- = 2J+ - 










(8.12) 



(8.13) 



we see that the results of (8.12) are all consistent with the correlators (5.12) predicted on 
general grounds. 



The remaining correlators with J are governed by the following expectation values. 
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As a result, we find the corresponding correlators, 
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Using again (8.13) and the vanishing of the correlators of mixed chirality, we derive, 

kb p+ 



{J+ip)J^ 

{r++ip)j+ 



-p)) 
-p)) 



A7TG3 p_ 

CyCE p^ 
AnGs p. 



(8.16) 
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These results hold to all orders in g+. The coefficient on the first line coincides with the one 
computed for the purely magnetic solution (since it is independent of g+), while the second 
coefficient is the one announced in (5.18) and (|5.19 ). 



8.2 Sourcing only T 



We source T by turning on g^^l while keeping g 
(7.9), the field Ni is turned on since s = s'^'L 



0. In view of 



(0) _ (0) _ (0) 

46(7^1. One might now solve successively 



the fluctuation equations of (6.30) and (6.31) following the scheme of (6.34). Actually there 
is a simpler and more illuminating way of proceeding using the exact SL(2, R) invariance of 
the reduced field equations in the far region, given in (2.9). In the far region, the source g^^l. 



may be turned on by performing the following SL(2, R) transformation A on the fields, 
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(8.17) 



Using the transformation rules of (2.9) on the full fields L, M, N, E, P,V of (6.27) leads to 



the following exact solution to the fluctuation equations in the far region, 

Mi(r) = Li(r) 
Ei(r) = Ni{r) 



/iMo(r) 








2/lLo(r) 



Pi{r) = -/l^o(r) 



lit 



This far region solution is incompatible with the matching conditions in the overlap region, 
however, because the non-zero source s[?L = 4:bg^^l_ ^ implies that t++ 7^ is non-local. 



in view of (7.11). Therefore it cannot be consistently omitted in our calculations. This 



shortcoming is easily remedied by adding to Mi the homogeneous solutions Lq and Lq, with 
coefficients to be determined by the matching in the overlap region. We will do so below. 

Since the above solution has a non-trivial Pi-field, we are inadvertently turning on an 
unwanted source in the field A_. From the preceding section, we know how to turn this 



source off by a shift in as in (8.3). Since the fluctuation equations are linear, the two 



solutions may simply be added in the far region. The resulting total solution takes the form. 
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(8.19) 
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while we continue to have Vi = 0. This solution holds to all orders in g+ as may readily be 
checked by explicit calculation. 



8.2.1 Absence of sources and matching the solution in the overlap region 

The absence of the sources gfl. and A^^\ ai 
following relations on the metric coefficients, 



The absence of the sources and A^^\ and the presence of the source g^}_ imposes the 
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and the relations 
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on the Maxwell coefficients. We solve (8.21) for r by exploiting the overlap region relation 
p+cr_ — p_cr+ = 0, then use this result for r to solve (8.20) for and we find, 
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The second relation provides 3++ with a simple expansion in powers of g+ given by, 

'+-f 

P- " 

n > 3 (8.23) 



.[2] 



sll = 

P- 



Using the remaining equations of (8.20), and the matching conditions for the metric (7.11), 

= ^^L°L t^ll = til = t^^- = (8.24) 
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Substitution of given by (8.23) into the expressions for immediately reveals that 
these quantities are local, and so of course is t^^_ . Retaining only source terms and non-local 
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contributions in the fields gives Vi = and, 
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(8.25) 



This solution holds to all orders in g+, within the approximations of section 5.3 Extracting 



the leading asymptotics, it is now verified that only the fluctuations Li and A^^i have a source, 
while Ml, and a± do not. 

The presence of a linear term in Li (r) as r — )■ oo implies that we have turned on a source 
to the operator T"*"" as well as to T . More precisely, we have introduced the coupling of 
g^^l. to the combination 

7 - 2(3r+- = 4r++ + 8/32r__ (8.26) 



From the sub-leading asymptotics in (8.25), we derive the non-local parts of the expectation 
values. We clearly find J_ 
of the following correlators, 
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T+_ = T = Tk = 0, and these relations imply the vanishing 
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These relations are consistent with the relations ( 5.10[ ) and (5.12) derived on general grounds. 
The remaining expectation values are non-vanishing, and given by, 
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whence we derive the expectation values of the normalized current and stress tensor. 
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These relations hold to all orders in g+. We readily deduce the correlators in momentum 
space, and find, 



cvCeP+ 
pI cWeP+ 



(8.30) 



Using the vanishing of the corresponding mixed chirality correlators from (5.10) and (5.12) 
we obtain our final result, 
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(8.31) 



The result on the first line agrees with the first line of (8.15): this gives a check that a 



correlator which may be evaluated in two different ways does indeed come out uniquely. The 
first term on the second line reproduces the purely magnetic result on the second line of 



(5.1), while the second term gives the result announced in (5.4) and (5.19). 
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9 Discussion 



The purpose of this work was to compute long-distance correlation functions in the charged 
magnetic background solution as a means of probing the low-energy dynamics of the dual 
system of fermions at finite density and magnetic field. The results obtained shed light on 
some properties of this theory, but various questions remain. 

In particular, the existence of a quantum critical point in this system is now well es- 
tablished, both numerically and analytically, but we lack a good understanding of what is 
driving this transition. Since we are working with a fully top-down construction in which 
the dual gauge theories are completely specified (e.g., A/'=4 SYM), we can in principle hope 
to obtain complementary descriptions of the transition mechanism within both gauge theory 
and gravity. It would be very instructive to have such an understanding of a finite density 
quantum phase transition from these two points of view. Note that in contrast to models of 
quantum criticality involving probe fermions [301 Ell 1321 EH] , here the gauge theory fermions 
are expected to be important players in the dynamics of the phase transition, rather than 
spectators. 

The results obtained here for the correlators of the stress tensor and U(l) current provide 
new information about the nature of the critical point. At least as probed by these operators, 
the low energy theory appears to retain the character of a 1+1 dimensional CFT, albeit 
with twisted Virasoro generators and renormalized propagation speed. This behavior seems 
consistent with a general 1+1 dimensional Luttinger liquid description. On the other hand, 
what is puzzling is that this gives no hint as to why the low temperature thermodynamics 
exhibits the nontrivial specific heat behavior displayed in equation (1.2) of the Introduction; 
such behavior is not characteristic of Luttinger liquid behavior, but rather indicates the 
breakdown of such a description. 

To appreciate this point, it is useful to compare the situation found here to one in which 
the low energy effective theory is Lorentz invariant. In a Lorentz and scale invariant theory 
with a traceless stress tensor, it is simple to show that there is a direct relation between 
the stress tensor correlator and the specific heat exponent; the power law obeyed by the 
former determines the latter. For instance, a 1/x^ falloff of the stress tensor correlator in 
1+1-dimensions fixes the specific heat to be linear in temperature. In our case, we also have 
a 1+1-dimensional effective theory with a stress tensor correlator exhibiting 1/x^ falloff, but 
the specific heat is not, in general, linear in the temperature. This can be traced back to 
the lack of Lorentz invariance; specifically, due to the fact that the stress tensor cannot be 
made to be both symmetric and traceless (see [M] for a recent discussion of this point). 
The conclusion is that in a non-Lorentz invariant theory one cannot necessarily expect to 
see nontrivial specific heat behavior mirroring itself in the behavior of the stress tensor 
correlator, which is consistent with our computations. 

The question thus remains as to what other observable might exhibit qualitative features 
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sensitive to the degrees of freedom that are going critical at the quantum phase transition. 
One natural possibility is to consider correlation functions of charged fields in the bulk. We 
expect that these can be computed using the same methods employed in this paper. Ideally, 
such correlators will exhibit behavior at the critical point that can be used to elucidate the 
mechanism driving the transition in the dual gauge theory. We leave these investigations for 
future work. 
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A Scaling laws for low temperature thermodynamics 



In this appendix we review the fc-dependent behavior of the low temperature thermodynamics 
around the quantum critical point. Recalling that we can, without loss of generality, assume 
A; > 0, there turn out to be three distinct regions : < A; < 1/2; 1/2 < A; < 3/4; and 
k > 3/4. 

We recall that for our system the quantum critical point arises at a particular value of 
the dimensionless magnetic field, B = Be, obeying the following properties. For B < Be the 
entropy density remains finite at zero temperature; for B > Be the entropy density goes to 
zero linearly with temperature; and for 13 = B^ the entropy density vanishes as a nontrivial 
power law, s ~ T°. The values of i?c and a are /c-dependent quantities in general. 

For A; < 1/2 there is no quantum critical point, as the zero temperature entropy density 
is found to be nonzero for any value of B. This was established numerically in [TU]. The 
special role of k = 1/2 can be seen by examining the near horizon solution (2.15). This 
solution clearly breaks down at k = 1/2, and for A; < 1/2 has the wrong large r behavior 
of g++ to match onto an asymptotically AdSs geometry. This near horizon geometry thus 
plays no role for A; < 1/2, and the low temperature behavior is instead governed by a finite 
entropy density solution whose form has not been determined analytically as yet. 

Next consider the range 1/2 < k < 3/4. The low temperature thermodynamics in this 



region is found to be controlled by the near horizon solution (2.15). This is a scale invariant 



geometry, and as discussed in section 4.1, the thermodynamic behavior associated with this 
scaling is: 



k 



a 



(A.i: 



In Table [T] we compare this prediction against numerical data, and the agreement is seen to 
be excellent. As we have noted, this scaling behavior only holds in the range 1/2 < k < 3/4, 



k 


a 


(1 - A:)/A; 


Aa 


0.55 


0.8182 


0.8182 


0.0003 


0.60 


0.6667 


0.6667 


0.0003 


0.65 


0.5389 


0.5385 


0.0003 


0.70 


0.4319 


0.4286 


0.0026 


0.749 


0.3606 


0.3351 


0.0032 



Table 1: Numerical results for the exponent a, its statistical standard deviation Aa, and 
comparison with the analytically predicted values (1 — A;)/A;, in the scaling relation s ~ T" 
for the range 1/2 < k < 3/4. 
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and we now discuss why this is the case. In order for a given near horizon, zero temperature 
solution, such as (2.15), to control the low temperature thermodynamics, it must be the 
case that there exists a finite temperature solution that shares the same asymptotics. Only 
in this case will the scaling behavior of the zero temperature solution manifest itself at 
finite temperature. What happens for k > 3/4 is that at finite temperature the function V 
appearing in the general metric Ansatz ( 2.5[ ) turns out to grow as one moves out from the 
horizon. This is incompatible with the asymptotics of (2.15) in which V is constant. The 
fact that this behavior sets in at A; = 3/4 can be seen numerically, but a simple analytical 
explanation is presently lacking. 

In the range A; > 3/4 the low temperature thermodynamics is not associated with the 
scaling behavior of a near horizon solution, but rather emerges from the more involved 
asymptotic matching analysis performed in [TT]. This analysis, which confirms the earlier 
numerical results presented in [10], leads to the scaling behavior s ~ T^/^. Note that this 
agrees with the value of a quoted in (A.l) for /c = 3/4. 



B Perturbation theory and WKB for scalar correlators 

In this appendix, we derive a systematic perturbation theory in powers of small ^ for the 
scalar correlators of section [3} and develop a WKB approximation for large ^. 

B.l Perturbative expansion in power of ^ 



For 1/2 < A; the term in equation (3.22) is suppressed compared to the term for 

large z. For small z, this situation is reversed, but the u"^ term then dominates over both 
and z'^~^'^, as long as < 1. Thus, in the interval 1/2 < < 1, one expects a convergent 
perturbative expansion in the z^~^^ term to hold. The starting point is equation (3.22), 



which we prepare in a form adapted to carrying out perturbation theory in powers of ^, 

xV" + V - (x^ + = (B.l) 

The solutions for ,^ = are the modified Bessel functions Iu{x) and Ky{x). Regularity of 
the ^ = solution y^o near the horizon as x — )■ cxd requires ^o{x) = Kv{x) up to an overall 
normalization. The Green function corresponding to the differential operator on the left 



hand side of (B.l) with these large x asymptotics is given by, 

G{x, y) = e{x - y)K,{x)h{y) + e{y - x)KMUx) (B.2) 
Using the Wronskian relation xK^I^ — xI^Kl = 1, it is readily checked that, 

{x^dl + xd^ -x^- G{x, y) = -y6{x - y) (B.3) 
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The differential equation (B.l), together with the horizon asymptotics may now be trans- 
formed into the following integral equation for y), 



^{x) = K,{x) -e I dy Gix, y) y'-'' ^{y) 



(B.4) 



Next, we show that the perturbative solution to (B.4) is point-wise convergent in the interval 
1/2 < A; < 1. The perturbative solution is obtained by iteration, 



n=0 



^n+i{x)=l dyG{x,y)y^-^''^M (B.5) 



where f^o^x) = K,y{x), and the Green function G was given in (B.2). To investigate the 
convergence properties of the expansion in powers of ^, it will be useful to have the asymptotic 
behavior of the Bessel functions. For a; — 0, we have, 



X 



.2j T{i^ + 1) 
while for X — )■ oo, we have, 



KJx) 



1 f X 



2 V2 



Tiv) 



(B.6) 



Iv(x) 



V2 



KJx) 



■KX 



y^2x/7r 



(B.7) 



Subleading terms will not be needed here. We begin by considering the first iterationp] 



^,{x) = K,{x) dyy'-^'Uy)K,{y) + Ux) dyy^-'^ K,{y f (B.8) 

JO Jx 

Clearly, for all 1/2 < /c < 1, both integrals are point-wise convergent for all finite x. 

To investigate higher order terms in the expansion (B.5 ), we need the asymptotic behavior 
of which may be obtained from the asymptotics in (B.6) and (B.7). For small the first 
integral in (B.8) behaves as x'*"*^'^, while the second integral behaves as follows. 



X 



-2u+4-4k 



if 



X 







if 



dyy'-''K,{y)' ~ 

Thus, as X — )■ 0, the behavior may be captured as follows, 

(pi{x) ~ x'K^{x) e = Min(2z/, 4 - Ak) 



2u + 4<4k 
2v + A>Ak 



(B.9) 



(B.IO) 



^^The case A: > 1 can be acommodated as well by relaxing the lower integration boundary to a non-zero 
positive value (such as oo), and absorbing an infinite term proportional to K^(x) into <^o(a^)- 
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which, for < 1, is softer than the leading order term Ky{x). For large x, the integrals 
behave as, 

/■oo 

/ dyy^-^'^K^ivY ~ x^-^'^e-^^ (B.ll) 

For k > 1/2, and as a; — )■ oo, the leading behavior arises entirely from the first integral. For 
k > 3/4, the asymptotic behavior is precisely that of the leading order term, 

'fi{x)r~.K,{x) (B.12) 

Thus, the perturbative expansion in this region of k is uniformly convergent. For k < 3/4, 
an extra power arises, and we have, 

^i{x) ^ x^-^''K^{x) (B.13) 

the behavior of which is harder than the lowest order K^{x) term. Still, the contribution 
to each order is point- wise convergent, and we easily derive the following general leading 
asymptotic formula for l/2</c<3/4, 

, / 3— 4fc \ 

This leading order asymptotic behavior may be resummed, and we find, 

2(3 -4A;) 

a formula which demonstrates that the perturbation theory is well-defined. 



ifix) ~ exp (- Sf 1 K{x) (B.15) 



B.2 WKB expansion 

A key observation throughout is that the powers of p± occurring in the expansion parameters 
are independent of u. This suggests that, for certain regimes, we should be able to take z/^ 



large and use WKB. To this end, we use equation (3.20) as our starting point since it clearly 
exhibits both free parameters p"^ and p^. We now make the following change of variables. 

In 2; = a; -|- In A (f){z) = (p{x) (B.16) 

where A is a constant. (The change of variables of ( B.16[ ) is not to be confused with the one 



of (3.21).) The scalar equation (3.20) is transformed into the standard Schrodinger form. 
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Although and appeared as two independent free parameters in (3.20), only the com 



bination ^ of (3.23) is intrinsic, while p merely sets a scale for the coordinate 2;, or shifts 



an origin for the coordinate x. The addition of the constant In A to x precisely allows for 
such a shift. We will choose A such that it sets the coefficients of both potential terms in 
(B.17) equal to one another, thereby shifting the natural origin to x = 0. The corresponding 



condition on A is as follows. 



A 



4fe-2 



36p+ 



(B.18) 



With this value of A, the common factor of the potential terms becomes. 
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1 



\2k 



(B.19) 



and the Schrodinger equation of (B.17) takes the form 



g2x ^ g(4-4fc) 



2 2 

V if 



(B.20) 



Since the objects we are after will be independent of z/, we take v large. Specifically, for 
E = —h'^u'^ 3> 1, this equation admits a perfectly fine WKB expansion, provided h <^ 1, 
or ^ ^ 1. This range of parameters includes, for example, the long distance limit in which 
p+ — )■ at fixed The WKB method involves a turning point at x = x.,,, given by. 



E 



(B.21) 



For X > X*, we have an exponentially decaying branch, while for x < x,,,, two oscillating 
exponentials contribute, whose behavior accounts for the z'^ and z~'^ asymptotic branches. 
The WKB method provides a systematic expansion in powers of h, independently of the 
value of E. Therefore, our results for the form of the expansion parameter should be valid 
for all E, and thus all u, including the physical region > 1 of our problem. 



C Linearized reduced equations in light-cone gauge 

The reduced fiuctuation equations in light-cone gauge around the full charged magnetic brane 
solutions were derived using Maple. They are relatively involved, and will not be presented 

here. In part, they may be solved iteratively. The component of the Einstein equation 

may be solved for hrr in terms of hx by, 

2 ~ 

hrr = — „ t2 ^K (C.l) 
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The r— component of the Einstein equation may be solved for hr+ in terms of by, 

hr. = --[-^h, + —h,j (C.2) 
The — component of the Maxwell equation may be solved for a+ in terms of dr by, 



a+ = -jy^ (^(KoLohr)' + 2khhr) (C.3) 



The H — component of the Einstein equation may be solved for in terms of hx and a. 
(the dependences of the other fields having been eliminated using the above relations) by, 



r 



-3LlKoKfhK + 24hK + SLliK'^YhK + QK^LoL'fhK 
-9KoK',LoL',hK + QKliL'^ffiK - ^iKlL,p,A!^h^ (C.4) 



This leaves only the r component of the Maxwell equation and the rr component of the 
Einstein equation for the fields hx and 5^. They satisfy the following pair of coupled second 
order equations, 

U - /l^ + —riK - -f^'^K - JJ^'^K + llK + -p^fiK 

Mopl- 2p+p_ ~ 8~ 2{K'^f ~^ , 8iK,A'^P- . 

H f^riK 73~"-K ~ 72"^ ^ ~^ ^7 

-^0 -^0 ^0 -f^o -J-^o 

1 _ _ 1 _ 1 

= --KoLQa'^ - 2KoLlL'Qa[. - -L^K'^a'^ - -K^LlL'^ar - kbLlL'^ar 

-K,Ll{L',ra^ + 2kh'^~a^ + ^^~a^ + ^-^~a^ - '-K,M,pl~a^ 

Kq ZKq Kq z 

+p+p_KoLoar - -LlL'^KQCLr - -LlK'^ar + ip-LoA'^hx (C.5) 
The remaining equations are automatically obeyed. 

D Gauge change in the overlap region: derivations 

In this Appendix, we present the derivation of the solutions in the overlap region to the 



gauge change equations of (6.23), which we repeat here for convenience, 

a^e^P- = Smc^p" + t/^FxM + (D.l) 
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Imposing rotation invariance of (6.3) on Kmn -, (^m ^ hMN, and om gives = v? = 0. We 
impose the conditions for light-cone gauge of (6.4) on Hmn and a^, and for LMN gauge of 
(6.6) on Hmni an d aM - We assign the expressions for the hght-cone gauge solution in the 
overlap region of (6.16) and (6.18) to the remaining components of Hmn and cim- Equations 
(D.l) for the remaining components of Hmn are, 



h++ 


(r) 


= 4i6p_| 


.ru (r) + 


8ibq / 2k 

— - ic+r^'^ + c_r 










(r) 


= 2bu' 


V) + 2ihp. 


Y.ru^{r) + 2ibp-ru 




(r) 


= Aibp- 








(r) 









)- 



Aq^r 



2^2fc-l 



u (r) — ^Tzr-, —u [r) 



2k -I ' k{2k-l 



k{2k-l) 



u r 



(D.2) 



while those for the gauge field are given by, 



a_(r) 



ip+^(r) + (2fc-l)^^^^^^^± gr'^-Vfr) 



= ip^U{r) 

The equations determining ^ and the remaining components of u^^ are given by. 



(D.3) 










■ 9 9i'— 1 

iq p_r 



u^'ir)' + ip^u^ir) + (r) — ^ . — - w [j-] — v+r" — f _r 

OK [jLirh J- j 



-l-o- 



ru [r) + 
ru'ir)' + 



zp„M (r j 
ip^u^ir) 



2i / \ 
(r)' (aw+r'^ - (a + l)w_r-^-'"j 



863r2 6A;(2A; - 1 

C+ fc-2 I ^^^-l-fe 

4 4 



^(r)' + ^r"-'' + —r-^-^ + qr''-^ u+{r) 



(D.4) 



From (D.2) and (D.3), it is clear that, in the overlap region, the metric and gauge fluctuation 
fields h and a are determined completely in terms of v±, c±, u'^,u^, and 6. In turn, vJ',u^ 
and 9 are governed by differential equations (D.4) which only involve the data f±, c±. 



D.l Form of the LMN gauge solution 

Since the overlap region is contained in the far region, it is appropriate to match the solutions 
in the overlap region to (D.2), (D.3), and (D.4) with the form of the LMN gauge solutions 
in the far region. The latter are already known from equation (4.3) of [11], and may be 
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parametrized as follows, 



h (r) = s r + t 



a+(r) = a^r+'' + a° 
a_{r) = r"'' + a° 



(D.5) 



where the coefficients w±, t+_, a^., and a]_ are independent of r, and the con- 

straint equation implies s+_ = 0. Terms with different functional dependence on r have 



been omitted from (D.5), since they may be neglected in the overlap region. 



From equations (D.2) and (D.3), we can infer which functional forms of m*" 



u 



and 6 are 



allowed to contribute to the functional forms of h±± , 
allowed functional forms are found to be as follows, 



hx and a± specified in (D.5). The 



e{r) 



^0 ^-1 ^-2k ^-2fc+l 



,1 „-2fc+l -2k+2 

5 ' 



(D.6) 



From these functional dependences, it is clear that all terms involving r'^ and r ^ " may 



r ± 
U , U 



be dropped from (D.4), since they can never contribute to the relevant functional forms of 



9 in the overlap region given in (D.6). 
Furthermore, by scaling a factor 



P+P- 



out of r, we see that several terms in (D.2), 

'^^^(P-)*^ compared to 

other terms in the same equation, and are thereby suppressed for ^ ^ 1. This is the case for 



(D.3), and (D.4) exhibit coefficients that contain powers of ~ (p 



the fourth term of the second equation in (D.2), and the fourth term in the ffist equation of 
(D.4), both of which are suppressed by compared to the second terms in each equation. 



It is also the case for the third term in the third equation of (D.4) which is suppressed by ^ 



compared to the second term, upon making use of the second equation of (D.4). 



Omitting all irrelevant terms leaves simplified equations for the metric fluctuations, 



Atop+ru (r j H c+r + c_r — 

p„ ^ ^ 2A; — 1 

26u^(r) + 2ibp^ru^ {r) + 2ibp-ru^ {r) 

4:ihp-ru^{r) 



2^2fc-l 



2^2A: 



u r — 



Aip^q r 

k{2k - r 



M"'"(r) 



(D.7) 



and for the gauge field fluctuations, 

= u^{ry + ip^u'^{r) + ip_u^ (r) 
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,, ip^u [r] 
= ru+ r '+ ^ ^ ' 



353^2 



= ru (r) + 



353^2 

= ^(r)' + ^r'=-2 + ^r-i-'= + gr'=-iM+(r) (D.8) 



Note that the fluctuation hx decouples from all the others. 



D.2 Solutions in LMN gauge for metric fluctuations 



To solve the simplified equations of (D.7) and (D.8), we begin by taking the derivative of 
the first equation in (D.8), and eliminating the derivatives of using the second and third 



equations. This gives the following equation for alone, 



= u'\r)" + 



453^3 



u''{r) 



(D.9) 



The second term may be neglected in the overlap region, since it is suppressed by a factor 



of p^/r ^ 1. We are left with ■u'" 



0. The solutions to this equation and to the second 



and third equations in (D.8) may be parametrized as follows. 



■u^(r) 
u~{r) 



Un + 



863r 
8b^r 



(D.IO) 



The integration constants independent of r. In the result for u an extra term 

proportional to Uq/t"^ has been omitted, since its functional dependence is not in the list of 
(D.6). Enforcing also the first equation of (D.8), and omitting a term suppressed by a factor 



of p /r, gives the following non-trivial relation between the integration constants. 



Ml 



(D.ll) 



We note that the functional dependences involving r which were listed in (D.6 ) for u'^, u^, 
do not actually occur. As a result, the last two terms of in (D.7) do not contribute. The 



term in c+ will also not contribute there because of its functional dependence, and neither 
will the term in c_ because it is momentum-suppressed in view of (6.17). As a result, the 



metric fluctuations are indeed given by (D.5) and (6.24), with the coefficients given in (6.25). 
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D.3 Solution in LMN gauge for gauge field fluctuations 

Using the above result for the solution for 9 in the overlap region may be readily com- 
puted, and we have, 



(D.12) 



where 9q is a further integration constant. Using this expression in the solutions for a± 
of (D.3), and retaining only the functional behaviors of the solutions (D.5) suitable in the 
overlap region, we find. 



.1 



6i6c+ (lP+{k — 1] 



a:, = {2k - 1) + ^^-f — + -ff — S++ 



a 



Ak 



(D.13) 



Within our low momentum approximation, we may further set c_ = 0. 
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